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INTRODUCTION  (by  J.D.  ZUND) 

The  major  part  of  this  research  report  is  essentially  the  Ph.D. 
dissertation  of  Wayne  Moore,  "Conformal  Geometry  and  Differential  Geodesy," 


§1 


Mr.*: 


useful  to  mathematical  geodesists  who  are  interested  in  further  exploring  the 
fascinating  ideas  of  Martin  Hotine.  In  addition  to  the  dissertation,  I  have 


effect,  it  sets  the  stage  for  the  topic  to  be  investigated  in  the 
dissertation.  The  chapter  concludes  with  biographical  material  on  Antonio 
Marussi  and  Martin  Hotine,  and  their  proposal  to  create  a  three-dimensional 
geodesy  in  the  discipline  which  we  call  differential  geodesy. 

Chapter  II  —  Conformal  Geometry  —  is  a  self-contained  introduction  and 
systematic  exposition  of  conformal  geometry,  i.e.,  the  geometry  of  a  pair  of 
N-dimensional  Riemannian  spaces  VN  and  VN  which  are  related  by  a  conformal 
transformation.  It  presumes  a  prior  acquaintance  with  tensor  analysis,  and 
can  be  regarded  as  supplement  to  the  discussion  in  HOTINE  [1969]  (Chapter  10, 
pages  55-62)  and  EISENHART  [1949]  (Chapter  II,  pages  89-95).  The  latter  is 
the  standard  reference  in  English  on  the  topic,  however,  it  is  very  incomplete 
and  really  merely  a  partial  summary  of  known  results.  Chapter  II  begins  with 
a  brief  excursion  into  the  modern  formulation  of  conformal  geometry  in  Section 
52.  This  is  admittedly  somewhat  abstract,  however,  it  is  necessary  to  make 
precise  the  notion  of  a  conformal  transformation,  i.e.,  mapping,  between  VN 
and  VN.  This  notion,  and  what  is  meant  by  adopting  the  same  local  coordinate 
system  on  both  VN  and  v  ,  is  obscurely  done  in  the  classical  literature, 
e.g.,  EISENHART  [1949].  Sections  §3-5  then  contain  a  systematic  development 
of  conformal  geometry  using  the  classical  formulation  with  particular 
attention  being  given  to  the  dimensions  N  =  2  and  N  =  3  which  are  of 
primary  importance  in  non-relativistic  differential  geodesy.  The  most 
interesting  mathematical  aspects  of  conformal  geometry  occur  when  N  >  4,  and 
the  cases  N  -  2  and  N  »  3  are  often  inadequately  treated  in  the 
literature.  Section  §3  introduces  the  basic  ideas,  concepts,  and 
tensor-theoretic  quantities  encountered  in  conformal  geometry,  and  Section  §4 
discusses  their  behavior  when  N  =  2  and  N  »  3.  Section  §5  considers 
integrability  conditions  and  conformally  flat  spaces.  The  presentation  of 
this  material  in  the  usual  references,  i.e.,  EISENHART  [1949]  and  SCHOUTEN 
[1954],  is  wholly  inadequate  and  riddled  with  errors.  Moreover,  the  original 


material  is  also  confusing.  The  key  results,  which  are  the  "X-Representation" 
of  the  curvature  tensor  (see  pages  29-30)  and  Schouten's  lemma  (see  Theorem 
5.3,  on  page  33),  were  originally  done  in  Schouten's  symbolic  notation  which 
was  not  widely  understood  and  which  Schouten  himself  later  abandoned.  The 
material  is  quite  subtle  and  intricate,  and  I  regard  Section  §5  as  being  the 
clearest  and  best  treatment  of  the  subject  which  I  have  seen.  The  topic  of 
conformally  flat  spaces  is  likely  to  be  of  considerable  importance  in  future 
work  in  differential  geodesy.  Both  HOTINE  (1966a,  1966b]  and  MARUSSI  [1985] 
(see  pages  169-176)  suggested  that  such  spaces  are  important  in  studying  the 
propagation  of  light  in  continuous  isotropic  refracting  media.  Virtually 
everything  in  Chapter  II  is  necessary  to  understand  the  conceptual  setting  for 
resolving  Hotine's  Conjecture. 

Chapter  III  —  Triply  Orthogonal  Systems  —  is  devoted  to  introducing  the 
basic  notions  occurring  in  Hotine's  Conjecture.  Section  §1  reviews  the  theory 
of  triply  orthogonal  systems  of  surfaces  and  presents  some  elementary  examples 
of  such  systems.  It  also  includes  the  Dupin  Theorem  (Theorem  1.1,  pages 
39-41)  and  the  Generalized  Dupin  Theorem  (Theorem  1.2,  pages  41-42).  The 
former  is  proven  in  Section  §1,  however,  the  latter  is  more  complicated,  and 
its  proof  is  deferred  until  Chapter  IV.  Section  §2  is  devoted  to  giving  a 
tensor-theoretic  derivation  of  the  Cayley-Darboux  equation  and  explicitly 
exhibits  the  general  form  of  this  equation.  The  material  in  this  section  has 
been  accepted  for  publication,  see  ZUND/MOORE  [1986],  and  will  appear  in  the 
near  future.  It  turns  out  that  the  Cayley-Darboux  equation  is  the  critical 
result  in  refuting  Hotine's  Conjecture.  As  will  be  shown  in  the  supplement, 
Hotine's  version  of  this  equation  is  wrong,  and  this  led  him  into  believing 
that  this  equation  was  always  identically  satisfied  in  a  flat  three- 
dimensional  Euclidean  space  E^ .  Section  l<3  contains  a  statement  and  proof  of 


i.L*, 


Liouvilie's  Theorem  (Theorem  3.1,  pages  46-51).  This  important  result  is 
discussed  (without  attaching  Liouvilie's  name  to  it)  in  HOTINE  (1969)  (see 
page  56  and  his  footnote,  and  also  in  HOTINE  [1966a,  1966b]).  The  theorem  is 
rather  deep  and  delineates  the  possible  types  of  conformal  transformations 
between  N-dimensional  Euclidean  spaces  E^  and  Ej^.  Despite  its  importance, 
this  result  is  rarely  proven  in  the  literature,  and  Section  §3  concludes  with 
a  self-contained  proof  which  is  based  on  that  given  in  BIANCHI  [1910].  A 
general  N-dimensional  argument  is  given  in  DUBROVIN/FOMENKO/NOVIKCJV  [1984], 
but  it  is  inconplete  and  is  not  changed  in  the  second,  1986,  Russian  edition 
of  this  book. 

Chapter  IV  —  Differential  Geodesy  —  presents  the  analytical  apparatus 
for  analyzing  the  Hotine  Conjecture  and  the  flaws  in  his  alleged  proof.  After 
formally  stating  his  conjecture  in  Section  SI  and  why  he  hoped  it  would  be 
true,  in  Section  52  the  question  of  isometric  immersion  of  a  surface  V ^  in  a 

is  discussed  together  with  the  behavior  of  under  a  conformal 
mapping.  Section  53  introduces  the  formalism  of  Ricci  rotation  coefficients 
and  indicates  how  they  are  affected  by  a  conformal  mapping  of  V^  into  V^. 
Although  in  HOTINE  [1969]  this  kind  of  reference  system  —  it  is  convenient  to 
call  it  a  triad  —  was  implicitly  employed,  this  important  topic  was  not 
explicitly  utilized.  This  was  a  serious  omission  in  Hotine's  analysis,  since, 
as  we  will  show  in  the  supplement,  it  will  allow  us  to  conclusively 
demonstrate  that  his  Cayley-Darboux  equation  is  wrong  and  not  an  identity!  We 
believe  that  this  triad  formalism  is  particularly  suited  to  the  requirements 
of  differential  geodesy  and  will  prove  to  be  an  important  part  of  future 
developments  in  mathematical  geodesy.  Section  54  uses  this  formalism  to 
concisely  establish  two  important  criteria:  the  Schouten-Eisenhart  Theorem 
(Theorem  4.1,  pages  59-61)  and  the  Ricci-Finzi  Theorem  (Theorem  4.2,  pages 
61-63),  for  the  conformal  flatness  of  a  V  .  Theorem  4.3,  (pages  63-64) 


furnishes  a  new  Schouten-Eisenhartlike  criterion  for  a  V^.  A  more  abstract 
version  of  the  results  of  Section  §4  has  been  submitted  to  Tensor  N.S.  for 
publication.  The  properties  of  congruences  of  curves  being  normal,  geodesic, 
and  canonical  are  studied  in  Section  $5  using  rotation  coefficients  and  a 
proof  of  the  Generalized  Dupin  Theorem  is  given  (pages  70-71).  This  new  proof 
has  been  accepted  for  publication,  see  MOORE/ZUND  [1986J.  An  example  of  how 
to  find  a  canonical  congruence  is  presented  in  the  Appendix  (pages  101-103). 
Section  §6  contains  a  rotation  coefficient  formulation  of  the  Cayley-Darboux 
equation.  Section  $7  now  begins  our  analysis  of  Hotine's  Conjecture  by 
outlining  the  steps  of  his  argument.  Section  58  shows  by  using  Theorem  8.1 
(page  81)  that  part  of  Hotine's  proof  is  false.  Section  59  examines  when  the 
conformal  image  of  a  normal  congruence  of  curves  can  be  a  geodesic  normal 
congruence  and  suggests  that  another  step  of  Hotine's  argument  is  shaky. 
Section  510  briefly  considers  the  question  of  a  conformal  mapping  between  a 
pair  of  surfaces  and  V.,,  and  shows  that  a  formula  for  the  conformal 
image  of  the  geodesic  curvature  given  in  MARUSSI  [1985]  (see  his  page  150)  is 
incorrect.  Finally,  in  Section  511  our  results  are  applied  to  give  a  critique 
of  Hotine's  argument.  It  is  shown  that  further  flaws  occur  in  his  procedure: 
he  specialized  the  conformal  function  defining  his  mapping  and  made  essential 
use  of  coordinates/equations  which  are  valid  only  at  a  point  of  a  V^.  The 
effect  of  either  of  these  mistakes  is  to  demand  that  his  (curved)  surfaces 
degenerate  into  (flat)  planes.  None  of  these  steps  can  be  easily  rectified, 
and  when,  taken  together,  they  strongly  suggest  that  Hotine's  argument  is 
fatally  flawed,  and  that  his  conjecture  is  false. 

The  Supplement  (written  by  J.D.  Zund)  now  furnishes  the  conclusive 
reasons  why  Hotine's  conjecture  is  false.  It  is  intended  as  a  companion  to 
our  joint  paper,  "Hotine's  Conjecture  in  Differential  Geodesy,"  which  will  be 
published  in  Bulletin  G<-od<-sique .  In  this  supplement  we  not  only  extend  the 


material  in  Dr.  Moore's  dissertation,  but  provide  the  details  of  our 
refutation  of  the  Hotine  Conjecture  which  had  to  be  omitted  in  our  paper. 
After  some  introductory  comments  in  Section  51,  we  translate  Hotine's 
Cayley-Darboux  equation  into  the  rotation  coefficient  formalism  and  prove  that 
it  is  not  equivalent  to  the  true  Cayley-Darboux  equation  in  Section  52.  This 
shows  that  the  Cayley-Darboux  equation  is  not  an  identity  as  Hotine  claimed, 
and  hence  his  conjecture  cannot  be  true.  Finally,  in  Section  53  we  discuss 
the  physical  consequences  of  this  result,  and  what  they  mean  for  differential 
geodesy. 

In  conclusion,  Dr.  Moore  and  I  would  like  to  express  our  gratitude  to  the 
Air  Force  Geophysics  Laboratory  and  in  particular  to  Dr.  Jekeli  for  his 
invitation  to  prepare  this  research  report.  The  cooperation  and  support  of 
this  research  under  contract  F  19628-86-K-0028  "Conformal  and  N'on-Conformal 
Transformations  in  Differential  Geodesy"  is  gratefully  acknowledged. 


Joseph  D.  Zund 
Principal  Investigator  and 
Professor  of  Mathematics 
New  Mexico  State  University 
Las  Cruces,  New  Mexico 
July  1987 
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CHAPTER  I 


INTRODUCTION 


hi.  GEODESY  AND  MATHEMATICS 

Theoretical  geodesy  is  the  study  of  the  size  and  shape  of  the 
Earth.  This  is  not  geographic  shape,  but  rather  the  shape  of  the 
surface  of  mean  sea  level  and  its  continuation  under  the  earth's 
crust.  This  surface  of  mean  sea  level  is  called  a  geoid  and  by 
definition  it  is  an  equipotential  surface  of  the  earth's 
gravitational  field. 

The  history  of  geodesy  may  be  divided  into  three  eras:  the 
spherical  era,  the  spheroidal  era,  and  the  geoidal  era.  The 
spherical  era  dates  from  the  time  of  the  Greeks  to  approximately 
1670.  The  spheroidal  era  stretches  from  1670  to  approximately  1830. 
The  geoidal  era  spans  from  1830  to  the  present.  An  eminent 
mathematician  was  pivotal  in  the  transition  between  each  of  these 
eras. 

Pythagoras  is  normally  credited  as  being  first  to  conceive  of 
the  earth  as  spherical.  He  reasoned  that  since  the  sun  and  moon  are 
spherical,  then  the  earth  must  have  a  similar  shape.  Aristotle  made 
what  was  probably  the  first  scientific  estimate  of  the  size  of  the 
earth.  However,  since  we  have  no  idea  of  the  size  of  his  unit  of 
length  (the  stadium)  it  is  impossible  to  estimate  the  accuracy  of  his 
work.  In  the  third  century  B.C.,  Erastosthenes,  considered  one  of 
the  founders  of  geodesy,  devised  the  idea  of  measuring  the  size  of 
the  great  circle  arc  between  the  North  and  South  poles.  His 


technique  was  based  on  measuring  the  difference  in  the  angle  of  the 
sun's  rays  at  Alexandria  and  Syrene  (now  Aswan),  which  he  assumed  to 
be  on  the  same  longitude.  Although  Erastosthenes'  estimate  is  15% 
too  large  when  compared  with  modern  estimates,  his  idea  is  very 
modern  in  spirit  and  the  error  is  due  to  the  crudeness  of  his 
measurements.  In  the  first  century  B.C.,  Posidonius  made  another 
estimate  of  the  size  of  the  earth  based  on  the  difference  in  the 
angle  the  star  Canopus  made  with  the  horizon  at  Alexandria  and 
Rhodes.  His  estimate  was  also  11%  too  large. 

After  Posidonius  interest  in  geodesy  lapsed  for  some  1500  years 
when  it  was  rekindled  by  the  need  for  accurate  maps  and  the 
exploration  of  the  New  World.  In  1617,  Snell  used  triangulation  as  a 
method  to  determine  distances.  This  was  a  breakthrough  since  this 
method  was  much  more  accurate  than  direct  measurement.  In  1669, 
Picard  used  a  telescope  to  determine  latitude  as  well  as  in 
triangulation.  Picard's  work  is  extremely  important  since  Newton 
used  his  estimate  of  length  of  the  arc  of  a  degree  of  latitude  to 
show  that  gravity  extends  beyond  the  surface  of  the  earth  and 
determines  the  motion  of  the  moon.  This  estimate  of  the  length  of  a 
degree  of  latitude  ultimately  became  the  basis  of  the  metric  system. 

The  period  between  Erastosthenes  and  Picard  is  the  spherical  era 
of  geodesy.  The  work  of  Isaac  Newton  inaugurates  the  spheroidal  era. 
Newton's  discoveries  in  mechanics,  i.e.,  laws  of  motion,  and  the 
formulation  of  the  law  of  gravitational  attraction,  are  crucial  steps 
in  determining  the  shape  of  the  earth.  Newton  used  a  theoretical 
argument  based  on  the  hydrostatic  equilibrium  of  the  oceans  to  show 


that  the  earth  was  an  oblate  sphere  with  the  major  axis  1/230  longer 
than  the  minor  axis.  The  modern  value  for  this  eccentricity  is 
approximately  1/290.  After  Newton  published  this  estimate,  many 
people  were  anxious  to  confirm  or  refute  it.  It  was  reasoned  that  if 
the  earth  were  flattened  at  the  poles,  a  degree  of  latitude  would  be 
shorter  near  the  poles  than  at  the  equator.  The  Cassinis,  a  family 
of  astronomers,  using  estimates  for  a  degree  of  latitude  in  the  north 
and  the  south  of  France,  determined  that  the  earth  was  a  prolate 
sphere  and  not  an  oblate  sphere.  This  announcement  naturally  aroused 
a  good  deal  of  controversy,  and  in  the  1730' s  the  French  Academy  of 
Sciences  sent  two  expeditions  —  one  to  Lapland  and  one  to  Peru  —  to 
settle  the  matter.  Their  measurements  showed  that  the  earth  was 
indeed  oblate  with  eccentricity  1/178.  It  is  interesting  that 
Newton,  sitting  in  his  room  in  Cambridge,  could  produce  a  better 
estimate  than  the  French  Academy  could  with  two  expeditions. 

With  the  initiation  of  mechanics  and  gravitation  by  Newton  other 
mathematicians  were  fast  to  extend  his  studies.  Euler  developed  the 
mechanics  of  rigid  bodies,  Lagrange  analytical  mechanics,  Legendre 
potential  theory,  Laplace  the  mechanics  of  rotating  fluid  masses.  A 
powerful  mathematician,  Clairaut,  computed  the  variation  of  gravity 
with  latitude.  Daniel  Bernoulli  and  Laplace  studied  tides  and 
methods  to  predict  them.  This  research  was  not  "pure"  research, 
which  was  then  applied  to  physical  problems,  but  mathematics  that  was 
developed  to  understand  physical  and  geodetic  phenomena. 

As  Newton  opened  the  spheroidal  era,  it  was  Gauss  (at  least  in 
spirit)  who  started  the  geoidal  era.  Gauss  became  geodetic 


consultant  for  the  Prussian  army  in  1799  and  in  1820  was  involved  in 
field  work  where  he  developed  new  instruments  for  surveying.  He 
adapted  his  method  of  least  squares  to  geodetic  measurements  and 
developed  the  Gaussian  probability  distribution  to  smooth 
observational  errors.  Using  data  from  geodetic  measurements,  he  was 
led  to  develop  his  theory  of  curved  surfaces.  Gauss's  intrinsic 
geometry  of  surfaces  eventually  led  Riemann  to  the  general  theory  of 
intrinsic  geometry. 

Mathematicians  and  physicists  who  followed  Gauss  developed  and 
extended  his  work.  Green  further  developed  potential  theory  (he 
coined  the  term) .  Stokes  calculated  the  undulations  of  the  geoid 
from  the  theoretical  ellipsoid  (Stokes  Theorem).  Rayleigh  and 
Poincare  extensively  studied  tides. 

Starting  around  1900  the  close  connection  between  geodesy  and 
mathematics  begins  to  diminish.  This  was  caused  by  exciting  new 
problems  in  mathematics  causing  it  to  become  divorced  from  physics. 
Moreover,  physics  was  evolving  into  modern  physics  and  leaving  behind 
classical  physics  with  which  goedesy  was  primarily  concerned. 
Moreover,  with  the  development  of  accurate  instrumentation  geodesy 
became  less  concerned  with  mathematical  methods  and  theoretical 
physics. 

^2.  HOTINE,  MARUSSI  AND  THREE-DIMENSIONAL  GEODESY 

The  technical  achievements  of  the  1950' s  made  it  possible  for 
geodesy  to  at  last  leave  the  surface  of  the  earth.  High  flying 
aircraft  and  satellites  could  measure  gravity  high  above  the  earth 


and  submarines  could  measure  it  deep  below  the  sea.  The  geodetic 
community  was  slow  to  take  advantage  of  these  new  technologies  as 
MARUSSI [ 1985) ,  p.  6,  notes: 

The  third  dimension  thus  appears,  in  practical  Geodesy,  as 
an  intruder  in  the  flourishing  paradise  of  the  Geodesist,  which 
comprises  the  surfaces  of  the  ellipsoid  and  the  geoid;  and  when 
this  intruder  comes  to  claim  his  proper  rights,  no  efforts  are 
spared  to  get  rid  of  him  as  quickly  as  possible,  and  with  the 
least  trouble,  by  means  of  a  weighty  battery  of  corrections  and 
reductions,  which  for  two  centuries  has  figured  in  every 
treatise  on  Geodesy. 

The  first  two  advocates  of  three-dimensional  geodesy  were 
Antonio  Marussi  and  Martin  Hotine.  Their  advocacy  was  for  a  point  of 
view  —  that  geodesy  is  inherently  three-dimensional .  Both  men  began 
to  apply  and  develop  mathematics  so  that  this  viewpoint  could  be 
practically  put  to  use. 

Antonio  Marussi  was  born  in  Trieste,  Italy  in  1908.  He  received 
a  Ph.D.  degree  in  mathematics  at  the  University  of  Bologna  and  then 
joined  the  Istituto  Geografico  Militare  (Italian  geodetic  and  mapping 
agency)  in  Florence.  During  his  twenty  years  with  them  he  modernized 
geodetic  procedures,  adopted  international  standards,  and  streamlined 
computing  schemes.  In  1952,  he  accepted  a  professorship  at  the 
University  of  Trieste. 

Marussi  introduced  advanced  mathematical  tools  such  as  tensor 
analysis  to  geodesy.  He  also  recognized  the  practical  aspects  of 
geodesy  by  doing  field  work  in  numerous  countries.  He  died  in  1984. 

Martin  Hotine  was  a  career  military  officer  and  rose  to  the  rank 
of  Brigadier.  He  attended  the  Royal  Military  Academy  at  Woolwich, 


England  and  was  commissioned  into  the  Royal  Engineers  in  1917. 
During  his  career  in  the  British  Army  he  pioneered  practical  methods 
of  topographic  mapping  using  aerial  photography,  did  geodetic 
surveying  in  East  Africa,  served  with  the  British  Ordinance  Survey, 
and  during  the  Second  World  War  he  served  as  Deputy  Director  of 
Survey  in  the  British  Expeditionary  Force.  Hotine  earned  a  degree  in 
engineering  from  Magdalene  College,  Cambridge,  during  his  time  in  the 
service. 

Hotine  retired  from  the  military  in  1946  and  became  the  first 
director  of  the  Directorate  of  Overseas  Surveys  where  he  served  until 
1963.  He  then  jointed  the  U.S.  Coast  and  Geodetic  Survey.  When  the 
Environmental  Sciences  Services  Administration  was  formed,  he  became 
a  research  staff  member  of  this  agency's  laboratory  at  Boulder, 
Colorado.  It  was  during  this  time  that  he  wrote  HOTINE[1969]  and 
pioneered  the  systematic  use  of  tensor  techniques  in  geodesy.  Hotine 
died  in  1968. 

Hotine  and  Marussi  came  to  geodesy  from  vastly  different 
directions.  Both  men  achieved  distinction  in  theoretical  and 
practical  geodesy  and  came  to  similar  conclusions  concerning 
theoretical  work.  Their  viewpoint  is  summed  up  by  MARUSSI ( 1985 1 : 

In  effect,  what  we  know  today  of  the  earth's  gravity  field 
owes  very  much  more  to  work  done  at  the  desk,  with  pencil  and 
paper,  than  to  observations  made  with  instruments  in  nature, 
(p.  5). 


§3.  NOTES  ON  THIS  DISSERTATION 

In  this  dissertation  we  have  examined  a  conjecture  by  Martin 
Hotine  on  triply  orthogonal  coordinate  systems.  In  Chapter  II,  we 
introduce  notation  and  develop  machinery  for  conformal  geometry.  In 
Chapter  III,  triply  orthogonal  coordinate  systems  are  discussed. 
Chapter  IV  then  deals  with  Hotine 's  conjecture. 

We  have  employed  a  somewhat  unusual  method  of  listing  items  in 
the  bibliography.  This  method  is  based  on  that  employed  by  J.A. 
Schouten  in  his  book  SCHO(JTEN[  1954  ] ,  which  is  probably  the  most 
comprehensive  book  on  tensor  analysis.  Items  in  the  bibliography  are 
listed  by  author's  name  and  listed  chronologically.  References  to 
specific  items  are  indicated  by  the  author's  name  in  capital  letters 
with  the  date  displayed  in  brackets.  In  case  of  several  items  in  a 
given  year,  the  dates  have  lower  case  Latin  indices  attached. 
General  references  to  an  author  are  given  by  citing  the  author's  name 
without  complete  capitalization.  We  are  grateful  to  Mary  Eberhardt 
of  the  Graduate  School  for  her  understanding  and  advice  on  how  best 
to  use  Schouten' s  system  in  our  bibliography. 

We  have  used  the  Einstein  summation  convention  on  repeated 
indices.  When  we  do  not  want  to  sum  on  repeated  indices  the  letters 
"NS"  for  "no  sum"  will  be  written  by  the  equation. 
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CHAPTER  II 


CONFORMAL  GEOMETRY 

*»1 .  INTRODUCTION  AND  BASIC  NOTATION 

The  purpose  of  this  Chapter  is  three-fold: 

(i)  to  introduce  the  notation  and  conventions  to  be  enployed, 

(ii)  to  relate  the  classical  and  modern  terminology, 

(iii)  to  provide  proofs  of  some  "well  known"  classical  results 
which  are  not  easily  accessible. 

Our  notation  is  essentially  that  of  EISENHART[ 1949] ,  with  slight 
modifications.  However,  our  presentation  is  considerably  more 
detailed.  Let  V.,  denote  an  N-dimensional  Riemannian  manifold. 
Initially  we  assume  that  N  >  3  and  the  metric  tensor  is  positive 
definite.  The  cases  N  =  2  and  N  =  3  will  be  dicussed  separately. 
The  symbolism  and  ”="  will  denote  "equal  by  definition"  and 

"identity",  respectively.  Free  indices  will  be  h,  i,  j,  k,  P,  m,  n 
and  summed  indices  (used  for  emphasis  when  only  some  indices  are 
summed)  will  be  p,  q,  r,  s,  t  .  In  the  case  N  -  2  or  3  no 
special  indices  will  be  used.  Later,  special  conventions  will  be 
made,  e.g.,  when  N  =  2  Greek  indices  will  be  used. 

Section  2  uses  its  own  notation,  which  will  be  introduced  there. 
Additional  notation  will  be  given  where  appropriate. 

‘<2.  THE  MODERN  FORMULATION 

Definition  (2.1).  Let  M  and  M  be  smooth  N-dimensional 
Riemannian  manifolds  (all  manifolds  are  assumed  connected  and  of 


class  C  ).  The  Riemannian  manifolds  M  and  M  are  locally 
conformal ,  wheneve  r , 

(i)  for  every  p(  N  and  p  €  M  there  exists  a  neighborhood 
of  p  :  U  C  M  and  a  neighborhood  of  p  :  U  C  M  ,  and  a 
map  f  :  U  ->  U  ,  f(U)=U  such  that  f  is  a 
diffeomorphism  with  f(p)  =  p  ; 

<n 

(n)  there  exists  a  C  function  X  :  U  -»  IR  such  that  for 
any  pair  of  tangent  vectors  Xp  and  Yp  in  the  tangent 
space  of  M  at  p  ,  T  (M)  ,  that  satisfies 

(a)  A  >  0  on  U  , 

(b)  for  the  inner  products  <  ,  >  and  <  ,  >_  in  the 

M  M 

respective  tangent  spaces  T  (M)  and  T. (M)  , 

P  p 

one  has 


(2.1) 


<f.xp.£*V:. '  MpI<VVm  • 


where  fx  :  T  (M)  -»  T.(M)  is  the  derivative  (or  Jacobian)  map  of 
*  P 

P 

f  .  The  map  f  is  a  local  conformal  map,  and  the  function  A  is 
called  the  conformal  factor. 

The  derivative  map,  fA  ,  is  defined  by 


(f*Xp)g  :=  Xp(g°f ) 


where  g  c  C  (M,(R )  . 


The  local  charts  { U , h }  and  (U,h)  of  M  and  M  , 


respectively,  are  related  by  the  commutative  diagram 
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and 


(2.5a) 


;  2 . 1 )  becomes 


<ei'V  :=  g(ei'ekJ 


gik 


(2.6) 


=  8  ^Ik  ‘ 


The  tensors  g  and  g  with  components  g  and  g  respectively, 
are  the  metric  tensors  of  M  and  M  ,  respectively.  This  is  the 
usual  tensor  expression  defining  a  locally  conformal  map  between  M 
and  M  (EISENHART! 1949] ) .  classically,  e.g.,  in  DUBROVIN/FOMENKO/ 
NOVIKOV] 1984 ] ,  one  often  encounters  the  expression 


(2.7] 


.  r  .  s 
dX  dx 

gi  ]  ~  grs  ,‘i  *  j 

J  tlx  dXJ 


which  reduces  to  (2.6)  by  virtue  of  our  definition  of  the  map  f#  in 
(2.4). 

Actually  (2.6)  is  improperly  written  since  it  relates  components 
of  tensors  defined  in  the  tensor  products  of  different  cotangent 
spaces.  The  correct  expression  requires  introducing  the  pullback  (or 
restriction)  mapping  of  the  cotangent  spaces,  i.e., 


f  :  T-  M  «  ...  «  T.(M)  -  T  (M)  »  ...  W  T  ( M) 

P  P 

P  P  K 


defined  by 


if  aHpHVj^,  ...,  vR)  =  a(  f  ( p) )  (  f  t  ( v^  ) ,  ...,  f  *  ( vk ) ) 
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where  a  is  a  covariant  tensor  of  order  k  and  v^  e  Tp(M)  ,  and 
requires  writing 


12.8) 


or  classically, 


(2.9) 


e*‘  2n 
f  g  =  e  g  , 


,*  2i> 

f  g.  .  =  e  g. 
’13 


ik 


However,  the  use  of  (2.6)  is  so  pervasive  in  the  literature  that  in 
this  dissertation  we  will  always  use  it.  We  will  always  use  the  term 
conformal  to  mean  locally  conformal  and  henceforth  write  M  -  VN  and 
M  -  VKl  as  in  EISENHART11949). 

We  now  give  expressions  for  the  Levi-Civita  connections  and  the 
curvature  operators  on  M  and  M  ,  respectively.  Writing  g  -  Xg  , 
,\  >  0  and  let  v  and  v  denote  the  Levi-Civita  connections 
compatible  with  g  and  g  ,  respectively.  Then  using 


(2.10)  <Vx,Y,Z>  *  j  (X<Y,Z>  +  Y<Z ,X>  -  Z<X,Y> 


+  <Z , ( X, Y] >  +  <Y,(Z,X]>  -  <X,[Y,Z]>)  , 


: 2 . 10)  is  just  the  definition  of  the  Christoffel  symbols)  we  have 


[2.11)  v^Y  =  vxY  +  |  ( (XfsjY  +  (Y+)X  -  <X,Y>v-k) 


where  41  :*  log  a  . 
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Calculating  the  curvature  form  for  v  we  have 


(2.12)  R(X, Y)Z  =  R(X,Y)Z  +  4  (h,(X,Z)Y  -  h,(Y,Z)X 
~~~  ~  ~  ~  Z  V  ~  ~  ~  •I'  ~  ~  ~ 


+  <X,Z>H  Y  -  <Y,Z>H  X)  +  i  [((Y*)(Z 
~  ~  'P~  ~  ~  'P-  4 

-  <Y,  Z>Mv^ll 2  )X  -  ((Xp)(Z4>)  -  <X,Z>  llv^ll2  )Y 


+  ((X*)<Y,Z>  -  (Y4/)<X,Z>v^)  ] 


where  R(X,Y)  is  the  curvature  form  of  v  ,  H  is  the  Hessian 

~  ~  'P  - 

tensor  which  is  a  tensor  of  type  (1,1)  on  M  with  H.X  :=  v  . 
-  'p~  x 

and  the  Hessian  form,  h,  ,  is  given  by  h(X,Y)  :=  <H,X,Y>  .  H  is 

\p  y  \p  \p 

self-adjoint  with  respect  to  the  Riemannian  metric.  See 
GROMOLL/KLINGENB ERG/MEYER [ 1968 ] . 

Our  viewpoint  and  analysis  are  always  local.  However,  the 
distinction  between  local  and  global  concepts  is  often  fuzzy  in  the 
geodetic  literature. 

A  stronger  notion  of  conformality  is  that  of  conformal 
equivalence . 


Definition  2.1.  Let  M  and  M'  be  N-dimensional  Riemannian 
manifolds.  Then  M  and  M'  are  conformally  equivalent  whenever 
there  exists  a  dif feomorphism 

f  :  M  -  M' 

such  that  f  defines  locally  conformal  maps  on  neighborhoods  of  M 


A 


and  M'  . 

2 

For  example,  consider  the  2-sphere  S  and  the  plane  E^  .  It 

2 

is  well-known  that  by  stereographic  projection  S  is  locally 

2 

conformal  to  ;  however,  since  S  is  compact  and  E^  is 

non-compact,  these  manifolds  cannot  be  dif feomorphic. 

5  3 .  THE  CLASSICAL  FORMULATION 

The  modern  definition  of  conformal  mapping  was  given  in  §2.  In 
sections  3-5  we  use  the  classical  definition  and  terminology. 

Definition  3.1.  If  the  metric  tensors  of  v.,  and  v  are  related 

-  N  N 

by 


3.1) 


2a 

g.  .  =  e  g.  . 


where  o  is  a  smooth  function  from  V^T  to  the  real  numbers,  then 
and  VN  are  said  to  be  conformally  related  or  just  conformal . 


NB:  The  choice  of  the  conformal  factor  as  e^7  instead  of  a 2 
or  ,\  (used  in  *i2)  is  mere  convenience.  The  choice  of  e20  is  nice 
for  differentiation  and  eliminates  many  unnecessary  factors  of  2 
and  ^  ,  but  has  no  geometric  significance.  The  function  o  in  e2*7 
is  the  conformal  function. 

We  now  define  the  Christoffel  symbols  of  : 

(3,2)  1  ijk  :=  2  <9ik|j  +  9jk|i  "  9ij|k)  ' 


16 


where  "  | "  followed  by  a  subscript  denotes  partial  differentiation 
with  respect  to  local  coordinates. 

By  direct  calculation  using  (3.1),  (3.2),  and  (3.3)  we  obtain 
the  Christoffel  symbols  for  VN  : 


(3.4) 


where 


:  3.5) 


.  k  .  k  j.  k  kp 

ii  ii  ii  yir  P 


and  where  we  have  written 


(3.6; 


a .  :=  a  i  . 

1  ll 


which  are  the  components  of  va  ,  i.e.,  the  gradient  of  a 


NB:  We  add  suffixes  omitting  a  differentiation  sign  only  on 

scalar  quantities  and  only  for  first  derivatives.  HOTINE[1969]  did 


not  follow  this  convention,  and  in  v  we  write 

N 


i  I  j  ='T|ill 


./  a.  .  . 

11 


We  will  denote  covariant  differentiation  with  respect  to  f\  ^  by 


and  with  respect  to  /  .  by 


.  Directional  derivatives 
will  be  denoted  by  "/"  ,  e.g.,  the  directional  derivative  of  f  in 


the  direction  A  is 


The  Beltrami  differential  parameter  of  the  first  kind  will  be 


denoted  by  A^a  where 

.  i  j 

A .0  :=  q  Ju.n.  . 

1  1  1 

The  Beltrami  differential  parameter  of  the  second  kind  (or 
Laplace-Beltrami  operator)  is  denoted  by  and  given  by 


(3.8) 


A,o  :=  gHj.  .  . 


2  i,j 

The  usual  Euclidean  expressions  for  A^o  and  A^o  are 


(3.9] 


AjO  =  (vo)*(vo)  , 


where  denotes  the  usual  Euclidean  inner  product  and 


(3.10) 


A^c  -  v* (vn )  =  v  o 


We  see  how  (3.9)  and  (3.10)  are  related  to  (3.7)  and  (3.8)  by 
denoting  the  Cartesian  metric  tensor  by 


(3.11: 

so  that 


g .  .  =('>.,  , 

yi3  13 


A  .0  =  ('> 1  .a  . 


w 

\-vy 

•Vj 

\s“S5 

elvSS 


Note  that  in  V. 


note  also  that 


=  gl3ai,j 


a .  .  ?  a  i  .  i  . 

i]  1 1 1  D 


a  .  .  *  n .  .  ; 

l]  1,3 


eara 

m 

AAS 


i 

V.v-.v 

>.V.J 


°lilj  =°il] ; 


and  by  definition  one  writes 


a .  ■  a  .  .  -  a  .a  . 

i]  1  3 


The  Riemann  tensor  of  will  be  denoted  by  >  the  Ricci 

tensor  by  R_  ,  and  the  scalar  curvature  by  R  .  The  Riemann  tensor 


of  VN  is  given  by 


(3.12) 


^hi jk  =  e  (Rhijk  +  ‘hijk' 


where 


(3.13) 


‘Shijk  ghkf'ij  +  gijnhk  9hjaik 
"  giknhj  "  Alaghijk  ' 


•  V  V  * 

i 

“•Tn-r 

%>.v 

•v  %*  %* 

•  vv 

fvvy 

>>v- 

<  */  v  * 

•  */  j» 

>w 

,vv 

VV'i 


H 


& 

v  v 

■vs*% 


%  v 
v  m*  V 

fttt 


>> 

.Vv 


9 


B 


3 


v 


i-' 


■s 


ghijk  ‘  9hj9ik  9hk9ij 


Some  useful  contractions  of  -hijk  and  9hij]<  are  9iven  bY 


(3.15)  S.  .  :=  Kh..  =  (N  -  2  )a  +  (A~u  +  (h  -  2 )  A  ^  )  g 
i]  i]n  13  ^  x  ■ 


or  alternately 


(3.16) 


2ij  =  <N  "  2)l°ij  +  gijAio]  +  2Tn=IT  Igij  ' 


where 


(3.17)  >;  :=  =  2(N  -  1 ) A +  (N  -  1)(N  -  2) 


(3.18) 


ID 

hk 


*  (1  -  N)9^  ' 


’hi  jk 


7i] 


(3.19) 


glDghkghijk  =  (1  "  N)N  ' 


(3.20) 


h  ,h,  .h 

9  ijk  j9ik  kgij 


6h 

D 


9ij  9ik  |  , 


(3.21) 


a  =  a 

ghijk  nigpj9qk 


9hj 

9hk 

9ij 

9ik 

and 


(3.22) 


hi  -  .hi  ,  h,  i  ,h,i 

g  =  :=  o^b,  -  o,.r>^  , 

jk 


i  ,  :=  <>  o,  -  11.  u  • 

3k  ]  k  k  3 
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and  where  h1  is  the  Kronecker  delta. 

3 

The  relations  ( 3.12)— { 3. 22 )  make  it  possible  to  calculate  the 

Ricci  tensor  and  scalar  curvature  of  V.,  .  The  Ricci  tensor  of  V  , 

N  N 

R^  ,  is  given  by 


(3.23)  R .  .  =  R.  .  +  (N  -  2)<>.  .  +  g.  .(A ,n  +  (N  -  2)A.ct)  , 
i]  l]  n  2  l 


1]  "1] 


and  R  ,  the  scalar  curvatures  of  ,  is  given  by 


(3.24)  R  =  e  (R  +  2(N  -  l)A,a  +  (N  -1 )  ( N  -  2)4^) 


As  usual 


(3.25) 


_  „h  hk_ 

Rij  ijh  ~  ^  ^ijk 


> 

3 


I 

>»J 

w; 

V  v\ 

■w  v 

\V  .A 

v/i 

'.i.V4 


(3.26) 


We  also  have 


R  :=  g1]R. .  =  R1. 

^  n  i 


(3.27)  g.  .R  =  g. .(R  +  2(N  -  l)Ayt  +  (N  -  1 ) (N  -  2)A.a)  , 
’ij  ’i]  2  1 


N  S  > 

)  _ 


from  which  it  follows  that 


(3.28)  g.  .R  -  g.  R  -  (N  -  1 ) (N  -  2)A,og  = 
^ij  i]  i  ij 

2  ( N  -  1 )  A  y’g  ■  • 

Z  1  J  / 


» 

Vv 


.  vv 

vv*. 

r.VV 


; 

►  v»»  « 

•yy, 

■>v-' 

v;o 

> 


.*  ».* 


.1  1.* 


hi]k  3hp  i]k 


^hijk  +  ( N— 2 )  ^ghjRik  ghkRij  +  ^ik^j  ^ij^hk' 


(N-1MN-2)  ghijk 


is  not  a  conformal  invariant  since 


(3.35) 


C  m  6  C 
''hijk  e  '“hi  jk 


The  tensor  is  called  the  Weyl  conformal  curvature  tensor. 
The  Weyl  tensor  is  completely  traceless,  i.e.,  all  contractions 
vanish.  It  satisfies  the  usual  symmetries: 


(3.36) 


(3.37) 


C*1  =  -C^ 

ijk  ikj 


Chijk  Cihjk  Cjkhi 


and  the  algebraic  Bianchi  identity. 


(3.38) 


ch.  +  Ch,  .  .  +  Ch ..  =0  . 

ijk  ki]  ;jki 


However,  it  does  not  satisfy  the  differential  Bianchi  identity.  To 
see  this  we  introduce  the  Finzi  tensor 


«»•»>  Lijk  Bij.k  -  *ik,j  *  '•«»*, j  -  »tjV  ■ 


which  satisfies 


(3.40) 


"ijk  “ikj  ’ 


23 


Then  the  conformal  analog  of  the  differential  Bianchi  identity  is 


(3.42) 


C  .  ..  „  +  C  .  „  .  .  + 
i]k,/’  if’}  ,k 


C  . 


j 


1  .  ,hr  ,hT  ,h,  th  rh 

. ,  ( o  .  L . .  .  +  n  L .  „  •  +  i>  „  L  . ,  +  q  . .  L  .  „  +  g .  „L,  ■  +  q .  . 
(N-2  j  lk/1  k  1  ^ j  /’  nk  ^ik  jP  ’iP  kj  nj 


The  singly  contracted  Bianchi  identity  is 


(3.43) 


R,  .  .  -  R.  ,  ;  + 


i j , k  ik,j  ^  mikj,h 

and  the  doubly  contracted  Bianchi  identity  is 


=  0 


(3.44) 


R1,  .  =  4  R  . 

k,  i  2  ,k 


From  this  it  follows  that 


(3.45)  L1^  =  0  . 

Moreover,  an  interesting  contraction  of  (3.42)  is  obtained 
(3.45) 


(3.46) 


ijk,h 


N-3 

N-2 


L 


ijk  ’ 


We  now  introduce  two  additional  tensors  —  the  Einste 
and  the  Cotton  tensor  —  which  are  important  tensors  in 
dimension  but  have  particular  significance  in  dimensions  2, 


by  using 


in  tensor 
arbitrary 
3,  and  4. 


The  Einstein  tensor  is  defined  by 


(3.47) 


S. 

1 


:=  R.  .  -  4  g 
3  13  2  y 


and  the  Cotton  tensor  by 


(3.48) 


A .  .  :=  R.  .  -  -j  q,  .R  . 
i]  4  ^i] 


**4.  SPECIAL  VALUES  OF  "N" 

We  first  consider  the  case  N  =  2  . 

The  elimination  procedure  used  in  §3  to  obtain  a  conformal 

curvature  tensor  breaks  down  for  R.  .  ,  R.  .  ,  R  ,  and  q.  R  when 

i]  i] 

N  =  2  .  We  first  show  that  every  V2  is  an  Einstein  space  (recall 
that  an  Einstein  space  is  a  Riemannian  space  such  that 


(4.1) 


R.  . 

i] 


R 

N  9 


*3 


or  in  the  case  N  =  2 


(4.2) 


R.  .  -  4  g.  .R  =  S.  .  =  0) 
13  2  ^13  13 


In  a  V 2  the  curvature  tensor  has  the  canonical  form 


(4.3) 


_h  . . , h  . h , 

R  .  =  K (t>  g  -  g.  .ft. 

ijk  ]^ik  ^13  k 


and  so 


(4.4; 


.  =  K( g .  .  -  2g .  . 


)  =  -Kg, 


i] 


Hence  we  have 


Therefore, 


(4.7) 


so  that  every  V.,  is  an  Einstein  space.  Now  observe  that  when 
N  =  2  we  have 


(4.8) 

R.  . 

=  R.  +  \-r>g-  •  , 

O 

1]  2  *i] 

(4.9) 

R  - 

e  ^  ’  ( R  +  2A  2«  )  , 

(4.10) 

g  ■  ■ 

=  g .  [  R  +  2A  ) 

’i]  2 

i] 

so  that  we  have 


(4.11) 


R  .  -  R.  -  4( g  R  -  q.  ,R)  =  0  , 

i]  i]  2  *i]  *i] 


or  using  the  definition  of  , 


(4.12)  S.  -  S  =  0  . 

i]  i] 

But  since  is  an  Einstein  space,  the  elimination  scheme  leading 

to  a  conformal  curvature  tensot  fails  at  this  stage.  This  is  why 
C^1  ..  and  C,  .  in  (3.33)  and  (3.34)  are  undefined  when  N  »  2  . 


On  the  other  hand  is  well  known  that  every  is  locally  conformal 

to  every  other  . 

Now  we  discuss  N  =  3  . 

The  dimension  N  =  3  is  much  more  complicated  than  the 
dimension  N  =  2  .  Nevertheless,  the  elimination  procedure  fails  to 
lead  to  a  conformal  curvature  tensor  and  in  fact  we  show  that  the 
conformal  curvature  tensor  discussed  in  *<3  vanishes  identically  in  a 


Before  we  show  that  the  elimination  method  fails  when  N  =  3  , 

we  derive  several  canonical  forms  for  the  curvature  tensor  in  a  V3  . 

The  simplest  derivation  of  a  canonical  form  for  in  a  is 

due  to  WILKES/ZUND[ 1978 ]  who  complete  a  problem  in  McCONNELL[ 1931 ] . 

i  i  k 

To  do  this  we  employ  the  Levi-Civita  dualizors  t  J  and  contract  on 
the  skew  symmetric  pairs  of  indices  in  the  curvature  tensor  Rhijk  » 
viz. 


(4.13) 


(ipq(  3rsR 


to  create  a  second-order  tensor. 


pqrs 

ipq,  jrs 


is  expressible  as  a 


3x3  determinant  in  the  metric  tensor  g^ 


(4.14)  fipq.jCS 


,P] 


,<33 


1  r 


1  s 


,ps 


Upon  multiplication  of  (4.13)  by  the  numerical  factor  of 
using  (4.14)  we  obtain 


1 

4 


and 


(4.15) 


S 


ipq,  DtsR 


pqrs 


where  S1-'  is  the  contravariant  Einstein  tensor.  McCONNELLf 1931 ] 
calls  this  tensor  the  Lame-  tensor.  However,  Lame's  work  was 
exclusively  three-dimensional  and  since  this  tensor  has  applications 
in  dimensions  N  ^  3  ,  the  name  "Einstein  tensor"  seems  more 
appropriate.  Equation  (4.15)  has  an  "inverse",  namely 


(4.16) 


R.  .  .. 
hi  jk 


'  hip  jkqR 


pq 


Equations  (4.15)— (4.16)  enable  us  to  write  down  several  canonical 
forms  for  Ry11jk  •  First, 


(4.17) 


Rhijk  ^hkSij  +  ^ijShk  ^hjSik 


c  _  K 

'  gik  nj  2  ’hi jk 


which  we  call  the  S-representation  °f  “hiik  • 
definition  of  the  Cotton  tensor  and  (4.18)  we  have 


Second,  from  the 


(4.18)  R, 


hi  jk 


g,  .  A  .  . 
’hk  l  j 


9-  A,  . 

’ij  hk 


^hj^ik 


^lk'^hj 


which  we  call  the  A -representation  of 
expansion  of  (4.18),  or  (4.19),  we  obtain 


R,  .  ..  . 

hi  jk 


Third,  upon 


R 

(4.19)  R,  .  =  q,  ,  R.  .  +  g.  R,,  -  g,  R.,  -  g.,R,  .  -  g.  .  ..  , 

hijk  ’hk  i]  ’i]  hk  ’hj  lk  ’ik  nj  2  ’hijk 

which  is  the  most  common  canonical  form  found  in  the  literature.  It 
shows  that  the  Weyl  tensor,  C^.  ^  ,  vanishes  identically  in  a  . 

However,  the  S-repres^ntati  n.  and  the  (-representation  are  more 


:  4 . 26 ) 


^hijk  e  ^  Rh  i  j  k  +  "'hi  j  k  ^ 


2n_  2tT 

=  e  R,  .  . ,  +  e  [q,  ,  (A  .  .  -  A .  . 
hijk  ^hk  i]  i] 


-  ^A.nq.  .)  -  g,  .(A.,  -  A., 

2  1  ^hj  lk  lk 


-  ^ikV*  -  gik(Ahj  -  Ahj 

-  phJAln]  -  e2nvghljk  • 


Therefore, 


(4.27)  R,  .  ..  -  g,  A.  .  -  g.  .A,.  +  g,  .A.,  +  g.,A,  .  = 
nijk  ^h]  i]  ^ij  hk  ^hj  lk  ^ik  hj 

e^r  (R,  •  ..  -  g,.A.  .  g.  .A..  +  g..A  . ,  +  q..A,  .) 

hi]k  ^hk  lj  yi]  hk  ^hj  lk  yik  hj 

but  this  reduces  to  0  =  0  by  the  A-representation  of  R^^ 
(4.18).  Thus  the  elimination  scheme  again  fails  to  give  a  conformal 
curvature  tensor  in  3  dimensions. 

The  foregoing  shows  not  only  that  the  Weyl  tensor  vanishes 
identically  in  a  ,  but  also  that  the  Cotton  tensor  and  the 

Einstein  tensor  naturally  occur.  The  former  plays  an  important  role 
when  we  consider  integrability  conditions  in  conformally  flat  spaces 
in  the  next  section. 


^5.  INTEGRABILITY  CONDITIONS  AND  CONFORMALLY  FLAT  SPACES 
Definition  5.1.  If  V_,  and  V,,  are  conformally  related  and  V„,  is 
flat,  i.e.,  is  locally  isometric  to  N-dimensional  Euclidean 

space,  ,  then  is  conform,!  1 1 y  flat.  A  conformally  flat  space 


if. 


4*. 


will  be  denoted  by  Cx,  . 

We  seek  to  develop  necessary  and  sufficient  conditions  for  a 
to  be  a  C  .  This  problem  will  be  addresed  in  two  parts,  first, 
N  >  3  ,  and,  second,  N  =  3  .  The  case  N  =  2  is  not  of  particular 
interest  since  as  noted  previously  any  is  conformal  to  any  other 

V2  • 

Preliminary  to  this  problem  we  suppose  N  >  2  and  assume  that 
the  curvature  tensor  of  VN  can  be  represented  in  the  form 

(5.1)  Rhijk  ~~  ^hkXij  +  ^ij^hk  ^hjXik  ^ikXhj 


where 

(5.2) 


and 

(5.3)  X  :=  gljX... 

are  to  be  determined.  The  representation  (5.1)  is  called  an 
X-representation  of  R^_k  •  From  (5.1)  we  have 


(5.4; 


R. 


13 


hk 

:=  g  R, 


'hi  jk 


=  NX.  . 
i] 


g.  .X 


-  X. 


13 


13 


so  that 


(5.5) 


R. .  =  (N  -  2)X. .  +  Xg.  .  , 
13  13  !3 


and 


(5.6) 


r  ;=  g^R  .  =  2(N  -  1  )X 


Hence,  we  obtain 


(5.7) 


1 


'ij  (N-2)  (Rij 


-  gijX) 


or 


(5-8)  Xij  N-2  (Rij  2 ( N— 1 )  9ij)  ' 

We  may  draw  the  following  conclusions: 

(i)  If  has  an  X-representation,  then  it  is  solvable 

for  the  X-tensor  and  its  contractions. 

(ii)  Expanding  the  X-representation  yields 

(5.9)  ^ijk  -  (N-2)  ^ghkRij  +  ^ij^k  ghjRik  ^ik^j 

_  R 

( N-l )  ghijk 

We  have  thus  proved 


Theorem  5.2.  For  N  >  3  the  Weyl  tensor  vanishes  if  and  only  if  the 
curvature  tensor  has  an  X-representation. 


Let  V„  and  v 
N  N 


be  conformally  related.  From  (3.12)  and 


(3.13)  we  have 


:  s.io) 


Rhi-ik  e 


2a. 


^iik  +  9hk°^  + 


hk 


%ini\ 


■  *  ».■  *2 

.V-. 


and  from  (3.23)  and  (3.24)  we  have 


(5.11)  R. .  -  R. .  +  (N  -  2 )a .  .  +  q.  .(A.o  +  (N  -  2)A.o 
lj  i]  n  ’it  2  1 


i]  "13 


(5.12)  R  =  e  (R  +  2 ( N  -  l)A2a  +  (N  -  1)(N  -  2)A]a) 


Hence  we  find  that 


(5.13)  g.  .R  -  g.  .(R  +  2(N  -  l)A,n  +  (N  -1 ) ( N  -  2)A,o) 
1]  1]  z  1 


Eliminating  A from  (5.11)  and  (5.13)  we  obtain 

(5-141  "ij  ‘  TiPT)  llRij  -  Rtj»  -  2T.PTT  -  V’1 

1 

-  g.  .A,o  = 

2  ^ij  1 

X.  .  -  X.  .  -  i  g.  .A.a  . 
i]  13  2  ^13  1 

Upon  replacing  this  value  of  in  (5.10)  we  see  that  R^^  *ias 

an  X-representation  if  and  only  if  R^^  ^as  an  ^-representation. 

Now  suppose  N  >  3  and  let  R^^  =  ®  •  Then  R^_.  =  0  and 
R  =  0  ,  so  that 


(5.15) 


a .  .  =  -X .  .  -  t  g. .A.o  , 
i]  i]  2  i 


(5.16) 


a .  .  =  a  .n  .  -  X.  .  -  -Jr  g .  .A  ,u 

1,3  13  13  2  ^13  1 


Replacement  of  this  in  (5.10)  shows  that 
has  an  X-representation.  That  is, 

Conversely,  suppose  that  c^ijk  =  8 
the  differential  equation  (5.16).  For 
symmetric  in  i  and  j  and  n 
integrability  conditions: 


Chijk  =  0  and  that  *hijk 
=0  in  a  CN  . 

and  that  a  is  defined  by 

a  to  exist  a.  .  must  be 
b] 

must  satisfy  the  following 


(5.17) 


i , j , k  i,k,  j 


=  a.  R .  .. 
h  i]k 


It  is  easily  seen  that 
(5.16)  we  obtain 


a 


i- 3 


is  symmetric,  and  differentiating 


(5'18)  °i,j,k  “  °i,k°j  +  V'j.k  '  Xij,k  -  2  ^ij^l^.k  • 
Interchanging  j  and  k  gives 

(5’19)  ai,k,j  =  °i,fk  +  °i"j.k  "  Xik,j  '  I  gik(Alo),j 
Thus ,  we  have 


(5.20) 


C  .  ,  —  CT  . 

1 ,  3  <  k  i,k,] 


ai,knj  -  °i,jak  +  Xik,j  "  Xij,k  +  2(gik(Alo),j  "  gij(Ala),k 


and  using  (5.16)  and  Theorem  5.2,  the  integrability  conditions  become 
merely 


5.21) 


X.  .  .  -  X..  .  =  0  . 

13- k  ik,3 


.  »v* 

v.v, 


v-y 

v  .sv 


\*  v". 

••  w 


» 


fejCi 


But  these  follow  from  (3.46)  when 


0  .  Thus,  a  function  a 


where 


satisfying  (5.16)  exists.  If  we  define  V.,  by  q.  .  =  e^°q.  . 

N  1  nj  ’ij 

<7j  is  the  function  of  (5.16),  then  a  straightforward  computation 
shows  that  =  0  •  Thus,  we  have  the  following  theorem. 

Theorem  5.3  (Schouten's  Lemma)  SCHOUTEN[  1921 1 .  A  V..(N  >3)  is  a 
CN  if  and  only  if  the  curvature  tensor  of  has  an 

X-representation . 

We  also  have 

Corollary  5.4.  A  ( N  >  3 )  isa  if  and  only  if  the  Weyl 

conformal  curvature  tensor  vanishes. 


In  S3  we  defined  the  Finzi  tensor: 


(5.22)  L.  ..  =  R.  .  .  -  R.,  .  + 


Jijk  *'i  j ,  k  l'ik,j  T  2 ( N- 1 )  (gikR,j  gijR,k‘ 
and  we  now  observe  that 


:  5.23) 


L.  ..  =  (N  -  2)  (X.  .  .  -  X.,  .  )  . 

i]k  i],k  lk,]' 


If  and  are  conformally  related,  then  the  Finzi  tensor  in 


V  is 

N 


(5.24) 


L .  ,  =  ( N  —  2 ) ( X .  .  ,  -  X. .  .) 

i]k  i] ;k  ik;] 


Writing  (5.24)  in  terms  of  un-hatted  objects  we  have 


=  L.  ..  +  (N  -  2)0  X. ,  -  a,  X.  .  +  Q . ,  g  a  X  . 

i]k  1  ]  lk  k  i]  ’ik’  p  ]m 


mp  . 

-q,  .q  vt  +  n  .  ,  -  a .  ,  .  -  o.o . 

i]  P  1’3>k  !.•«,]  k  l] 


+  °fik]  +  (^ik^P< 


d  a  .  - 

P  3m 


which  reduces  to 


9ij9  %%!>>  +  V’(9ikaj  -  9i fy) 
+  I  ^ij(4la),k  -  I  9ik(Alff,,j1  ' 


(5.26)  L.  =  L.  ..  +  O.X.,  -  a.  X.  ■  +  g.,gmpf'  x. 

i]k  i]k  ]  ik  k  lj  ’ik’  p  jm 


-  q.  .q  V!  X,  +  o  .  .  .  -  a .  .  .  . 

P  km  1,3, k  i»k,] 


By  the  integrability  conditions  (5.17)  on  n  we  obtain 


(5.27)  L.  =  L...  +  a .X. ,  -  n.  X.  .  +  g. .  gmprj  X . 

ljk  i]k  i  lk  k  i]  ’ik’  p  3m 

up  „  „h 

-  g.  g  ro  X.  +  H  ...  . 

’ij’  p  km  h  13k 

Hence,  if  has  an  X-representation,  then  L.^  =  Lijk  '  anc* 

if  Lijk  =  Lijk  '  ('*ien  Rhijk  ^as  an  X-representation.  Thus,  we 
have  proven 


Theorem  5.5  (Finzi's  Theorem)  FINZI [ 1922 ] .  In  V^'N  ■  3)  the  Finzi 
tensor  is  conformal  invariant  if  and  only  if  the  curvature  tensor  has 
an  X-representation. 


V'Vts.t'V,."*'  .*■ ,".y» 


Corollary  5.6.  In  V^(N  >  3)  the  Finzi  tensor  is  a  zero  tensor  if 

and  only  if  v._  is  a  Cv,  . 

1  N  N 


We  now  consider  the  case  N  =  3  . 

By  Theorem  5.5  we  know  that  the  Finzi  tensor  is  a  conformal 
invariant  in  since  the  curvature  tensor  of  always  has  an 

X-representation,  viz  ,  the  Cotton  tensor.  If  we  define 

a  by  (5.16),  then  the  integrability  conditions  are 

(5.28)  L .  =  0  . 

i]k 

If  we  insert  (5.15)  into  (5.10)  we  have 


(5.29,  Bhilk  -  0  • 

If  R,  .  .,  «  0  ,  we  have  R.  .  *  0  and  R  =  0  and  so  X.  .  ■  0  . 
ni  jk  ij 

Therefore,  we  find  that 


(5.30) 


=  0  . 


Hence,  we  have 

(5.31) 


since  L 


i  jk 


Theorem  5.6. 


is  conformal 


A  V3  is 


invariant.  Thus  we  obtain 


conformally  flat  if  and  only  if  the  Finzi 


tensor  is  a  zero  tensor. 


NB:  There  is  no  Weyl  tensor  in  a  . 

To  summarize,  we  have  obtained  necessary  and  sufficient 
conditions  for  a  v^(N  2  3)  to  be  conformally  flat.  Additionally, 
we  have  derived  an  often-neglected  result  (Theorem  5.5)  on  the  Finzi 
tensor . 


CHAPTER  III 


TRIPLY  ORTHOGONAL  SYSTEMS 


U.  TRIPLY  ORTHOGONAL  SYSTEMS  OF  SURFACES 

Since  a  point  on  a  surface  is  determined  by  two  parameters,  a 
point  in  space  may  be  determined  by  three  parameters  or  curvilinear 

coordinates.  In  E^  ,  for  example,  we  may  transform  our  Cartesian 

i  ,123.,  ,  ,, 

x  =  (x  ,x  ,x  )  by  means  of  the  equations 


(1.1)  x1  =  x1(u,v,w)  i  =  1,  2,  3 


where  u1  =  (u,v,w)  are  curvilinear  coordinates.  The  surfaces 


u  =*  ,  v  =  ,  w  =  C^  ,  where  C^,  C^,  are  constants,  are  the 
coordinate  surfaces  and  they  intersect  in  coordinate  lines.  Examples 
are  cylindrical  and  spherical  coordinates,  i.e., 


M  1  2.3 

(1.2)  x  =  u  cos  v  ,  x  =  u  sin  v  ,  x  =  w 


,  1  2  .3 

(1.3;  x  =  u  cos  v  cos  w  ,  x  =  u  cos  v  sin  w,  x  =  u  sin  v  , 


respectively.  The  line  element  in  our  space  will  then  take  the  form 


(1.4) 


,2  dx  dx-^  ,  m ,  n 

ds  =  (> .  .  — -  -  du  du 

ij  ,  m  ..  n 
J  d u  du 


Of  particular  interest  are  those  coordinate  systems  for  which 


(1.5) 


,  dX  dx^  - 

h .  . - =  0  n  /  m 

in  .  m  ,  n 
J  du  du 


h  *  I  ' 


t. 


tfhen  u  <  a  (1.9)  gives  ellipsoids,  when  a  <  u  <  b  it 

2  2 

represents  hyperboloids  of  one  sheet  and  when  b  <  u  <  c  it 
represents  hyperboloids  of  two  sheets.  If  we  re-write  (1.9)  as 


(1.10) 


2  2  2 
x  Y  z  _  . 

2  +  ,2  +  2  1 

a  -u  b  -u  c  -u 


2  2  2 
x  Y  z  _  . 

2  +  ,  2  +  2  1 

a  -v  b  -v  c  -v 


2  2  2 
x  y  z  . 

~ 2  +  2  +  "1  1 

a  -w  b  -w  c  -w 


2  2  2 
u  <  a^  <  b^  <  c 


2  ,.2,2 
a  <  v  <  b  <  c 


2  L2  2 
a  <  b  <  w  <  c  , 


the  line  element  takes  the  form 


.2  .2.2  ,  2 .  2  .  2 ,  2 

ds  =  h^  du  +  1^  dv  +  h^  dw  , 


where 


(1.11) 


,  2  _  1  .  _ ( u-v)  (u-w) 

“1  _  4  '  ?  2  2 

(a  -u)(b  -u)(c  -u) 


.  2  _  1  .  _ ( w-v)  ( u-v) 

n?  _  4  2  2  2 

(a  -v)(b  -v)(c  -v) 


,2  1  .  ( u-w )  ( v-w ) 

n3  "  4  2  2  2 

J  a-wl(b-w)lc-w) 


We  now  give  a  proof  of  one  of  the  fundamental  theorems  on  triply 


orthogonal  systems. 


Theorem  1 . 1  (Dupin):  The  curves  of  intersection  of  a  triply 
orthogonal  system  of  surfaces  are  lines  of  curvature  on  each  of  the 
surfaces.  See  McCONNELLf 1931 ] ,  p.  216. 


Proof.  For  clarity  we  present  the  proof  in  the  form  of  four 
assertions. 

Assertion  ( i ) :  We  have 

(1-121  ’  V"/  ‘  T  *  I  “■'>  *  l-2 

where  t  is  the  qeodesic  torsion,  t  is  the  torsion,  h  „  is  the 

g  3  up 

Euler  tensor,  0  is  the  angle  between  the  unit  surface  normal  and  the 
unit  principal  normal  of  the  curve  and  Xn  are  the  surface 
contravariant  components  of  the  unit  tangent  vector  to  the  curve. 

For  a  proof  of  this  assertion  see  McCONNELL[ 1931 ] ,  p.  214. 
Assertion  (ii):  If  a  curve  is  the  intersection  of  two  surfaces 
which  cut  at  a  constant  angle,  then  the  geodesic  torsions  of  the 
curve  on  the  two  surfaces  have  the  same  value. 

Proof:  Since  the  unit  surface  normals  along  these  curves  make  a 
constant  angle  the  rate  of  change  of  the  angles  0  ,  0  with  respect 
to  arclength  the  unit  surface  normals  and  the  unit  principal  normal 
must  be  the  same ,  hence 


Theorem  1 . 2  (Generalized  Dupin):  A  necessary  and  sufficient 
condition  that  a  third  family  of  surfaces  can  be  associated 
orthogonally  —  viz  be  a  component  of  a  triply  orthogonal  system  — 
to  a  given  pair  of  orthogonal  families  is  that  every  surface  of  the 
two  given  families  intersect  each  other  in  a  line  of  curvature. 

We  postpone  the  proof  of  this  theorem  until  Chapter  IV  S5,  but 


we  note  that  Generalized  Dupin  implies  Dupin  and  that  Dupin  furnishes 
the  necessity  of  Generalized  Dupin. 


S2.  THE  CAYLEY-DARBOUX  EQUATION 

Let  E^  denote  a  three-dimensional  Euclidean  manifold  having 
the  line  element 


ds^  =  q.  .du^du^  , 

yi] 


where  q.  .  =  h.h.  are  the  curvilinear  scale  factors  and 
’13  1  3 

u1  =  iu,v,w)  are  local  curvilinear  coordinates  which  are  functions 

^  12  3 

of  the  Cartesian  coordinates  x1  =  (x  ,x  ,x  )  .  Following 

DARBOUX[ 1910 ]  the  equations  u  =  ,  v  =  ,  w  =  ,  where 

C1#  C3  are  constants,  are  said  to  define  a  Lame  system  of 

surfaces,  i.e.,  a  triply  orthogonal  system  ,  in  and  each 

equation  defines  a  Lame  family  of  >  . 

It  will  be  convenient  to  work  with  the  functions  u,v,w  :  E^  — *  IR 

and  to  denote  partial  differentiation  with  respect  to  the  coordinates 

x1  by  subscripts,  viz.  u.  ,  u.  .  =  u..  ,  etc.  We  will  use  Cartesian 
j  e  1  13  31 

tensor  notation:  Latin  subscripts  range  from  1  to  3  and  repeated 
subscripts  obey  the  Einstein  summation  convention.  The  Schouten 
bracketing  convention  is  employed  to  denote  symmetrization  and  skew- 
symmetrization  on  subscripts. 


A  version  of  the  following  material  will  appear  in  Tensor,  N.S.  See 
ZUND/MOOREJ 1986 ] . 


/■AV 


The  orthogonality  conditions  for  the  curvilinear  coordinated 


surfaces  can  then  be  written 


(2. 2)  g12  =  Uiv.  =  0  ,  g23  =  v.w.  =  0  ,  g^  =  u.w.  =  0  . 

It  is  important  to  observe  that  the  conditions  in  (2.2)  are  symmetric 
under  a  cyclic  permutation  of  the  parameter  functions  u,  v,  and  w  . 
We  extend  this  observation  to  the  symmetry  principle:  for  any  partial 
differential  equation  obtained  by  dif ferentiation  of  (2.2)  with 
respect  to  there  are  two  analogous  equations  which  arise  by  a 
cyclic  permutation  of  u,  v,  and  w  . 

Suppose  that  the  Lame  family  is  defined  by  the  surfaces 
u  =  .  Choose  any  one  of  the  equations  (2.2)  involving  u  ,  e.g., 
uivi  “  0  '  and  differentiate  it  with  respect  to  x-3  : 


[2.3) 


u.v. .  +  u. .v.  =  0  . 

ii]  i]  l 


We  note  that  (2.3)  allows  us  to  replace  the  second  partials  of  v 
with  second  partials  of  u  and  we  call  this  the  replacement 
property. 

Contracting  (2.3)  with  w_  and  employing  the  symmetry  property, 
we  obtain  two  analogous  equations.  Upon  adding  two  of  these  and 
subtracting  the  third  we  have 


(2.4)  v..w..u..  +  u,.w.,v..  +  U..V..W.  =  0 

(i  j)  i]  [i  ]]  i]  [13)13 


But  since  this  expression  is  symmetric  in  the  subscripts  i  and  j  , 


it  reduces  to 


Now  calculating  the  partial  derivative  of  (2.5)  with  respect  to  x  , 
contracting  the  result  with  u^  ,  and  using  the  replacement  property 
on  the  partial  derivatives  of  v  and  w  yields 


(2.6) 


A.  .v .  w  .  =  0  , 
i]  i  ] 


where 


A..  :=  u  u. .  -  2u.  u. 

ij  s  ijs  is  js 


Equation  (2.6)  may  be  regarded  as  a  linear  equation  in  the  six 


unknowns : 


(2.8) 


Xij  "  V(iWj> 


and  to  emphasize  this,  it  is  convenient  to  re-write  (2.6)  in  the  form 


(2.9) 


A.  .X.  .  =  0  . 
1]  O 


If  the  unknowns  X^_.  can  be  eliminated,  we  will  obviously  obtain  a 
third-order  partial  differential  equation  involving  the  parameter 
u  .  To  effect  this  elimination  by  a  determinant  we  require  five 
additional  equations  involving  the  X^  .  It  turns  out  that  these 
are  easy  to  get.  First,  by  using  (2.8),  (2.6)  can  be  rewritten  as 


(2.10) 


u.  .X.  .  =  0  , 
i] 


and  by  using  (2.8)  we  also  have 


(2.11) 


X.  .  =  0  . 
i  i 


The  remaining  three  equations  aio  given  by 


by  virtue  of  the  orthogonality  conditions  (2.2).  Thus  the  desired 
elimination  is  given  by  the  determinant 


(2.13) 


A11 

A22 

A33 

2A23 

*31 

*12 

U11 

U22 

U33 

2u23 

2u31 

2u12 

1 

1 

1 

0 

0 

0 

U1 

0 

0 

0 

u3 

u2 

0 

u2 

0 

u3 

0 

ul 

0 

0 

u3 

u2 

U1 

0 

>>> 


Sks 

>>> 

£v>5> 

MV* 

>.v.^ 

"  V* 

\  \ 
v  „*  * 

v\'\' 

;■.->> 

-/••vv: 


which  is  the  Cayley-Da rboux  equation.  Substituting  (2.7)  we  may 
write  (2.13)  as  a  difference  of  determinants 


(2.14)  0  =  u 


,  v  «•  / 


ulls 

u22s 

u33s 

2u23s 

2u31s 

2u12s 

\  ;.v! 
.*•  *  V- 

U11 

u22 

u33 

2u23 

2u31 

2u12 

yv> 

1 

1 

1 

0 

0 

0 

• 

t.'v'  . 

U1 

0 

0 

0 

U3 

U2 

i  *  *  ■  m  * 

v'/-. 

v  *.  /. . 

0 

U2 

0 

U3 

0 

“1 

»  - 

0 

0 

U3 

U2 

U1 

0 

ulsuls 

u2su2s 

u3sU3s 

2u2sU3s 

2u3s“ls 

2u1sU2s 

U11 

u22 

U33 

2u23 

2u31 

2u12 

>&<< 

• 

1 

1 

1 

0 

0 

0 

U1 

0 

0 

0 

U3 

U2 

y&: 

0 

U2 

0 

U3 

0 

U1 

1 

0 

0 

u3 

u2 

U1 

0 

N-.o 

v\ 

•*>> 


It  should  be  noted  that  these  equations  show  that  the  Cayley-Darboux 
equation  is  linear  in  the  third  partial  derivatives,  cubic  in  the 
second  partial  derivatives,  and  quartic  in  the  first  partial 
derivatives  of  the  function  u  .  An  explicit  form  for  the 
Cayley-Darboux  equation  can  be  exhibited  by  expansion  of  the 
determinants  in  (2.13)  or  (2.14)  and  the  results  may  be  written  as 


where  the  summation  sign  denotes  the  sum  of  a  cyclic  permutation  of 
the  Greek  indices  n,  p,  ->  over  1,  2,  3.  Inspection  of  (2.15) 
shows  that  the  Cayley-Darboux  equation  contains  324  terms.  The  most 
general  example  given  in  DARBOUX[ 1910 )  contains  only  four  terms. 


h  3 .  THE  THEOREM  OF  LIOUVILLE 

Using  the  tools  we  have  developed  in  this  Chapter  and  in  Chapter 
II  we  give  a  concise  proof  of  one  of  the  most  amazing  theorems  in 
geometry  and  analysis. 


Theorem  3 . 1  (Liouville).  The  only  conformal  maps  of  E^  to 
E^(  N  >  3)  are  similarily  transformations  (isometries  and  homothetic 
maps)  and  transformations  by  reciprocal  radii  (inversions  in  a 
sphere) . 


This  is  in  sharp  contrast  to  the  case  N  =  2  where  there  is  a 
rich  supply  of  conformal  maps.  We  give  a  proof,  due  to  BIANCHI [ 1910 ] 
using  the  Lame  equations,  that  is  valid  in  three  dimensions  only. 
For  a  proof  for  general  N  see  DUBROVIN/FOMENKO/NOVIKOV( 1984 ] . 


Proof:  In  a  triply  orthogonal  system,  u1  =  (u,v,w)  in  E^  the 
line  element  has  the  form 


where 


2  2  2  2  2  2  2 
ds  =  h^du  +  h2  dv  +  h3  dvr 


ds  =  h,du 
u  1 

ds  =  h-dv 
v  2 

ds  =  h..dw 
w  3 


are  the  respective  arclengths  along  the  coordinate  lines. 

The  Lame  equations  are  S1^  =  0  and  following  McCONNELLf 1931 ] , 
p.  156,  we  have 


where  h  ,  :=  h  and  h  ,  : =  h  . ,  ,  .  We  also  have 

ab  a  I  b  abc  a  |  b  [  c 


(h.  h2l’l  +  (h,  h12 1 2  +  2  h13h23  "  0 

1  2  h^ 

(3.4)  (^-  h32)2  +  (p-  h23)3  +  —2  h21h31  -  0 

2  3 

(h7  h13 } 3  +  (h7  h31)l  +  ~2  h32h12  =  0 

3  1  h2 


-Vv'i 

'.v.v.'i 


•*.  <v 

e‘j  V  4 
I  < 


The  Lame  equations  are  necessary  and  sufficient  that  the  quadratic 


form  (3.1)  be  reducible  to 


(3.5) 


2  1  2  2  2  3  2 

ds  =  (dxxr  +  (dx  )  +  (dx  ) 


v,  i  ,12  3. 

where  x  =  (x  ,x  ,x  )  . 

Let  x1  =  ( x1 , x2 , x3 )  be  a  local  coordinate  system  in  E3  ,  with 
x1  =  (x3,x^,x3)  the  coordinates  of  the  image  of  x3  under  a 


*.■••.**  ■ 

f-  r<'' 
•V 

► 

■'>. 


conformal  map,  viz. 


(3.6) 


‘i  A.  1  2  3. 

X  =  f  (x  ,x  ,x  )  . 


rrr: 
%  % 
\  «\ 


Then  the  conformality  of  the  mapping  E^  — •  E^  requires  that 


.  A 

<>v  v' 


ratio 


■  --'.v  ■ 


\A> 


a\'  irj  T\nrjmj  1TV 


*J  PJIMR'W  WWWW I 


j*uw»v*WTiFwinu*wiwi 


WWWKTWWWWWjl 

?sd 


(3.7) 


(ft .  _.dx^dx^  J/iA^dx^dx^  ) 


i] 


be  independent  of  dx1  and  dx1  ,  viz. 


(3.8) 


ft.  .dx^dx^  =  (I/A^ja.  .dx^dx^ 
n  13 


where  A  is  a  function  of  x3  ,  A  =  A(x\x^,x3) 


In  the  Lamr  equations  (3.3)  and  (3.4),  putting  x1  =  u1  and 


each  h.  =  l/A  ,  we  obtain 


A23  0  ' 


(3.9) 


A31  "  0  ' 


A 12  0  ' 


where  A  ,  :  =  A  ,  ,  and 

ab  |  a  |  b 


(3.10)  An  +  A22  -  A 22  +  A33  -  A33  +  A 13  -  AAr\ 


where  A^  is  the  first  Beltrami  differential  parameter  with  respect 


to  g .  .  =  A .  , 

13  13 


Equation  (3.9)  implies 


(3.11) 


A  =  XU1 )  +  Y(x2)  +  Z(x3) 


and  substitution  into  (3.10)  gives 


•2  -2  -2 

(3.12)  X  =  Y  =  Z  =  rr  :  v.--  =  k  =  constant 


2  ( X  ^  y  +  z) 


where  "•"  denotes  differentiation  with  respect  to  the  independent 
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■  ',C  /  «'  .*  »'***•  ,*■„»*  ,*  •  "  •  ,*  ~ m-  '  •  *#-  '/  *  «  '4»  *  _  -  '  «  '  »  ’  .  '  »  "  *  *  •  *  a 


imi 


s^v-v. 


■  V  J.*.  .* 


/-  .>  .-a  1 


.... , 

^  %"  V 


*CvC vlv 


•-  --V  .> 

'•  .N 


.x  .'• 


\  A  . 


• _ _ 

>  \*s>  vl 

V*  v  V  Jl 


>\d 


« 


*vV:^ 


•  •  J 


/  «•  J 


variable. 

Two  cases  occur  according  to  whether 
(i)  k=0  or  (ii)  k^O. 

In  (i),  k  =  0  implies  X  =  Y  =  Z  =  0  hence  X,  Y,  and 
constant,  and  therefore 


(3.13) 


X  =  const. 


This  is  the  case  of  a  similarity  transformation,  and  f 
is  an  isometry  or  homothety. 

In  (ii)  it  is  convenient  to  write  k  =  2/c  where  c 
constant.  Then  integration  of  (3.12)  yields 


v  1  ,,  1  1  >2  .1. 

X  =  —  ( ( x  -  a)  +  b  ) 


(3.14) 


v  1  ,,  2  2.2  .2, 

Y  =  —  ( ( x  -  a)  +  b  ) 


Z  =  i  ( (x3  -  a3)2  +  b3) 


But  (3.12)  also  requires 


(3.15)  (x1  -  a1)2  +  (x2  -  a2)2  +  (x3  -  a3)2  +  b1  +  1 

3  .  .  , 

'  (x‘  -  a1 r  , 


(3.16) 


b  =  0 


Thus, 


<■.  /.  A  A  -■  v  v.  -■  .-.V  v. 


Hence , 


....  2  .  . 

(3.19)  6  .cb^dx3  =  ~  A i  .dx1^3 

so 

(3.20)  x1  =  — 

r 

which  are  the  formulas  for  transformation  by  reciprocal  radii,  i.e 
inversions.  Hence  the  theorem  is  proved. 
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CHAPTER  IV 

DIFFERENTIAL  GEODESY 


hi.  INTRODUCTION 

Martin  Hotine  in  HOTINEf 1966b J  and  in  HOTINE(1969]  sets  forth 
the  hypothesis  that  any  sufficiently  smooth  function  <p  :  ->  IP 

with  non-vanishing  gradient  can  be  a  member  of  a  triply  orthogonal 
coordinate  system  for  E^  .  Hotine  notes  that  this  is  at  variance 
with  classical  differential  geometry,  which  states  that  for  a 
function  to  be  a  member  of  a  triply  orthogonal  system  it  must  satisfy 
the  Cayley-Darboux  equation  of  Chapter  III,  '<2.  Hotine  then  asserts 
that  this  equation  is  an  identity  when  the  Lame  equations  are 
satisfied.  See  Chapter  III,  '<  3 .  Although  it  is  difficult  to 
determine  when  the  Cayley-Darboux  equation  is  satisfied  we  shall 
develop  machinery,  in  particular,  Ricci  rotation  coefficients,  that 
will  help  in  this  question. 

Hotine  wants  his  assertion  to  be  true  since  in  geodesy  the 
primary  object  of  study  is  the  shape  of  the  geoid.  If  its  potential 
function  were  a  member  of  a  triply  orthogonal  system,  it  would 
automatically  provide  a  natural  physical  coordmatization  of  E^  for 
the  geoid.  Unfortunately,  as  we  will  show  in  this  Chapter,  Hotine' s 
argument  is  false  and  his  hope  for  such  a  result  cannot  be  realized. 

k<  2  .  CONFORMAL  MAPPING  AND  ISOMETRIC  IMMERSION 

Throughout  this  Chapter,  Greek  indices  range  from  1  to  2  and 
Roman  indices  from  1  to  3 . 

Let  be  a  3-dimensicnal  Rremannian  space  with  metric  tensor. 


g^  ,  curvature 
symbols  I\_.  . 
tensor 


tensor  R,  .  ,  Ricci  tensor  R.  .  , 

ni]k  i] 

Let  be  conformally  related  to 


and  Christoffel 
with  metric 


(2.1: 


2o 

g.  .  =  e  g.  . 


where  a  :  -•  (R  is  a  smooth  function.  Let  be  a  2-dimensional 

subspace  of  with  metric,  or  first  fundamental  tensor. 


(2.2) 


i  j 

a  ,  =  g.  .x  x.J.  , 
aft  n] a  ft 


where  x^(u^,u^)  is  a  parameterization  of  ,  and 


(2.3) 


i  _  Ox 

a  a 

flu 


Equation  (2.2)  requires  that  V2  be  isometrically  immersed  in  . 
The  obvious  question  is  whether  under  the  conformal  map,  is 
isometrically  immersed  in  ? 

First,  we  have 


(2.4) 


aft 


i  l 

=  g.  x  x.J. 
i]  «»  ft 


so  that  and  V 2  are  conformally  related.  N.B.  :  As  in  Chapter 

II  in  defining  the  conformal  map  we  impose  the  same  local  coordinates 
on  both  v.j  and  ,  i.e.,  x1  =  x1  . 

Let  v2  have  unit  normal  vector,  f  ,  with  covariant  components, 

£ ^  ,  and  contravariant  components,  f1  .  Under  the  conformal  map  we 
have 
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By  Gauss's  formulae  we  will  deduce  the  second  fundamental  tensor  b 
of  V2  .  Using  the  notation  of  Chapter  II  to  denote  covariant 
differentiation,  we  have 


(2.7) 


rn,  n  n  cn 

P  b  ,  =  x  =  x  -  ?  a  a  , 
o  ,/<  /f  o/i 


(2.8) 


^  /  n  —  ri  v 

f  b  =  x  .  -  £  a  a  .  )  , 
o/i  (7  ,/i  s  /f  nji 


and 


(2.9)  s\p  =  fn(b^  -  a/s  a,,,)  =  e°f n(b(/j  -  u^a^)  , 


where  a  =  n  .F 
/f  i 


Hence,  we  obtain 


[2.10) 


bnh  ~  e  (bafi  n/£an(1) 


We 


must  now  check  that  the  formula  of  Weingarten  holds  in 


l  .e. , 


(2.11) 


cn  /»■>.  n 

f  .  =  -a  b  x 

;o  o/i  i 


By  direct  calculation  we  have 


so  that 


where  on  :=  g^o .  , 
y  1 


;  2 . 13) 


rn  -o,,.n  n. 

f  =  e  f  +  a  /cx 
;a  ,<>  /f  a 


However,  we  see  that 


(2.14; 


r  n  -a  fh  , ,  .  n 

-a  b  „x  =  -e  a  b  „  -  o  a  ,)x 
n[i  t  af)  /l  op  y 


-o  n 

-e  (f  +- 

,o 


n. 

ox) 

/f  « 


-o  ,  n  n . 

e  f  +  o  x  ) 
,o  /f  o 


and  thus,  Weingarten's  equations  are  satisfied. 

We  now  check  that  the  equations  of  Gauss  and  Mainardi-Codazzi 
hold  in  Vj  ,  i .e. , 


(2.15) 


=  b  bn 


"crfi'id  “ay^Pfi  “(ti  “nfi  +  ^Vlijk 


b„  b 


h  i  j  k 
x  x„xJx,  , 

O  I  i  "T  o 


(2.16) 


,  „  rh  i  i  k  n 

b  „  -b  „  -  R  .  . .  6  x  x,,x  =  0 
aP;y  ay ;[)  nijk  a  p  y 


Equation  (2.15)  is  the  Gauss  equation  and  (2.16)  are  the  equations  of 
Mainardi-Codazzi.  We  verify  these  in  the  usual  way  using  (2.2), 
(2.10),  (2.11),  and 


12.17)  *»  „  -  x"  -  R"  ,,  -  a"\nR  „ 

<r;/?p  op;/*  o/n  o  rurfh 


Employing  (2.10)  and  (2.11),  (2.17)  becomes 


letter  denoting  the  vector.  The  letters  a-g  will  denote  ennuple 
indices  and  h-p  will  denote  tensor  indices. 


Definition  3.1.  The  Ricci  rotation  coefficients  are  given  by 


1  i 

■>  ,  :=  e  .  ,e,  eJ 

abc  ai,j  b  c 


NB:  The  -i  ,  are  invariants  under  coordinate  transformations  but 
abc 

not  under  ennuple  transformations. 


The  rotation  coefficients  are  interpreted  as  the  rate  of 
rotation  of  the  ennuple  vector  e^  with  respect  to  e^  and  ec  . 

The  rotation  coefficients  are  related  to,  but  not  identical  to  the 
Cartan  structural  coefficients. 

It  is  easy  to  see  that 


(3.4) 


7  ,  =  0  , 

aab 


:#:j 

m 


(3.5) 


7  ,  +7  =  0  . 

abc  bac 


'.V.  -V« 

## 

>  V  V  V 

p 

,V 

*V  *0  -v’ 

i 

<>v. 
>„>■%> 
>.'V 
.*  v  v'V 
AV- 

>>>. 


,v.vv 

<0 


rf-.v.v; 


Under  a  conformal  map  V,  -*  V,  the  components  of  {e  } 

5  i  ~cl 


transform  as 


,v  A 

..V'.T-VV 


(3.6) 


e  .  =  e  e  .  , 
ai  ai 

i  -a  i 
e  =  e  e 
a  a 


.AJv 

O  V 
*•  V  , 


where 


is  the  same  function  as  in  (2  1).  We  now  compute  the 
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rotation  coefficients,  i  ,  of  the  ennuple  { e_^ }  in 


direct  calculation  we  have 


l  * 

e.  .=e...-/..e. 
ai,j  ai|]  lj  ak 


,  a  .  n,  k 

=  ee  . )  |  .  -el.,  e  . 
ai  ]  i]  ak 


a,  ..  k  ,k 

=  e  (o  .e  .  +  e  .  i  .  -  i  •  •  e  ,  -  o  .a  .e  , 

3  ai  ai  1 3  1  j  ak  1  3  ak 

tk  km  , 

-  6j"ieak  +  9ij9  Vak1 

a ,  m. 

=  e  (a  e  .  +  e  .  .  -  a  .e  .  -  a .  e  .  +  g.  .n  e^  ) 

j  ai  ai,3  3  ai  1  aj  ’13  m  a 


=  e(e.  .  -  o.e.+q.o.), 
ai,3  1  a3  ^13  /a7 


where  a  .  :=  o.e1  .  Thus,  we  have 

/a  1  a 


,7)  e.  .  =  e°  ( e  .  .  -  o.e.  +  q.  .0.) 

ai ; 3  ai,3  1  a]  *13  /a 


Moreover,  we  have  the  equation 


(3.8) 


1  -a,  1  1  .1  . 

e  .  =  e  (e  .  -  o  e  .+  o  .0  ,  ) 
a  ,-3  a  ,3  a3  3  /a 


by  virtue  of 


(3.9) 


gi];k  ■  0  ■ 


(3.10) 


91J;k  =  ° 


Hence,  we  find  that 


(3.11) 


abc 


=  e  .  .e, 
ai;j  b 


J  =  e~ 


(■» 


abc 


-  o  ,6 

/b  a  c 


0.6.) 
/a  be 


and 


(3.12) 


i  ,  =  e  i  .  , 

abc  abc 


when  a,  b,  c  are  distinct. 

We  now  give  a  definition  that  we  will  need  in  the  rest  of  the 
Chapter. 


Definition  3.2.  A  congruence  of  curves  in  is  a  family  of  curves 

such  that  one  of  the  curves  of  the  family  passes  through  each  point 
of  a  chart  of  V.,  . 

If  f  is  a  vector  field  defined  on  chart  of  ,  then  the 

integral  curves  of  f  define  a  congruence  of  curves  and  the  vectors 
f  are  the  tangent  vectors  to  the  curves  of  the  congruence. 

h  4 .  OTHER  CRITERIA  FOR  CONFORMAL  FLATNESS1 * 

In  Chapter  II,  $5,  we  showed  that  for  N  >  3  ,  VN  is  a  C^  if 

and  only  if  the  Weyl  tensor  C^^  vanishes;  and  that  for  N  >  3 

is  a  C..  if  and  only  if  the  Finzi  tensor  L.  vanishes.  In  this 
N  1  J  K 


1A  version  of  the  material  in  this  section  has  been  submitted  to 
Tensor,  N.S. 
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'.y.: 


.  A  AV.V 


section  we  will  give  criteria  for  conformal  flatness  involving 
orthogonal  ennuples. 


Theorem  4.1  (Schouten-Eisenhart)  SCHOUTEN[ 1921 ]  and  EISENHART[ 1949 ) . 
is  a  C  (N  >  3)  if  and  only  if  for  every  orthogonal  ennuple 


h  i  j  k  n 
.  e  e,  e  Je  .  =  0 


*abcd  ‘  Rhijkea  eb  ec  ed 


for  distinct  values  of  a,  b,  c,  d  =  1 ,  2,  . . . ,  N 


Proof.  The  first  half  is  easy.  If  v„  is  a  C.,  then 

1  N  N 


chijk  -  0  so  thlt 


_  _  h  i  i  k  „ 

C  ,  j  :=  C,  .  ,e  e,  e  Je,  =  0  , 
abed  hi]k  abed 


and  we  have 


(4.3)  Rabcd  N-2^acRbd  ^ad^bc  +  9bdRac 

“  gbcRad’  “  (N-2)(N-l)(gacgbd  "  9ad9bc)  ' 

and  for  an  orthogonal  ennuple  we  have  g  ^  =  6  ^  .  Thus  we  have  that 

R  ,  ,  =  0  for  distinct  values  of  a,  b,  c,  d  . 
abed 

To  establish  the  converse  we  must  first  exhibit  the  ennuple 
components  of  the  Weyl  tensor.  For  distinct  values  of  a,  b,  c,  d 
these  are 


Cabcd  Rabcd  ' 


\  V  /  * 


V/  w. 


(4.5) 


(4.6) 


Cabad  Rabad  +  N-2  Rbd  ' 


Cabab  _  Rabab  +  N-2  ^^b  +  Raa^ 


/->.■ 

w; 

Xttv 

&£* 

>.W 

!v.y> 

v>> 

vV/, 

fe's' 


(N-l ) (N-2)  ' 


where 


(4.7) 


R  ,  :=  R.  .e  xe,  3 

ab  lj  a  b 


Hence,  we  have  the  equations 


(4.8) 


(4.9) 


R  ,  ,  =  0  , 

abed 


*abad  N-2  Rbd  ' 


(4.10)  R  ,  ,  =  -  (R.  ,  +  R  )  +  vT/k, 

abab  N-2  bb  aa  (N-l) (N-2) 


These  are  easily  seen  to  be  equivalent  to 


(4.11) 


(4.12) 


(4.13) 


R  ,  ,  =  0  , 

abed 


Rabad  Rcbcd  ' 


Rabab  +  Rcdcd  Racac  +  Rbdbd 


from  which  it  follows  that  all  the  components  of  C  ^  ^  are  zero. 

We  now  give  a  more  geometric  criterion  for  to  be  a  C’N  . 
This  was  first  established  by  RICCI [ 1918 )  for  N  =  3  and  generalized 
by  FINZI [ 1922 )  for  N  >  3  . 


y.v;j 

•VyN 

cy.v:. 


I 

• . *.  i 

••/**.*•*  *i 
•  -  * 


V-  ■ 


.v'.v£ 

•  h  *  _  ■  w 
•  -  ■ 


v  VV’ 

This  was  first  established  by  RICCI  [  1918 )  for  N  =  3  and  generalized 

/Vv  1 

by  FINZI (1922)  for  N  >  3  .  ’ttSV 

£vi>^ 
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►  I 

$££ 

*.*.*••  V 


r-f  rr  rx  vx  yw\tw  vrw  wwwi ^rw.  v  V  V.  viniVTin  w  V 


Theorem  4.2  ( Ricci-Finzi ) .  is  a  ,  (N  /  3)  ,  if  and  only  if 


there  exists  an  orthogonal  ennuple  of  { e  J  such  that  the 

~  9 


correspondence  congruences  of  curves  (/ 


curves  having  tangent 


vectors  e_ )  are  normal  and  isotropic,  i.e.,  their  rotation 

~  9 


coefficients  i  ^  satisfy  the  respective  conditions: 


(4.14) 


•>  ,  =  0  , 
abc 


(4.15) 


abb  acc 


for  distinct  values  of  a,  b,  c  =  1 ,  2,  . . . ,  N 


The  condition  (4.14)  for  normality  imposes  N(N  -  1)(N  -  2)/2 


conditions  on  the  i^c  and  has  the  usual  meaning,  viz.  that  there 


exists  an  N-tuply  orthogonal  system  of  ( hyper ) surfaces  {>)  } 

9 


i-hogonal  to  the  respective  congruences  {/ 


The  isotropy 


seditions  '4.15'  impose  N(N  -  2'  additional  conditions  on  the 


and  by  the  familiar  geometric  interpretation  of  the  rotation 


■fficients  this  means  that  each  I  is  turning,  i.e.,  "rotating", 

9 


■he  same  rate  :n  each  of  the  directions  (b  ?  a)  specified  by 


:  he  the:  ronut  uences  I.  .  Thus  (4.14)  and  (4.15)  reduce  the  number 

b 


:f  independent  rotation  coefficients  from  N  (N  -  l)/2  to  N  . 


Ric^i  proved  this  for  N  =  3  by  directly  integrating  the  system 


if  equations  given  by  '4.14)  and  (4.15)  and  noting  that  the  resulting 


ipace  was  a  C,  .  For  N  >  3  ,  Finzi  observed  that  (4.14)  and 


(4.15*,  when  substituted  in  the  lotation  coefficient  expressions  for 


abed  ' 


v.v, 


£  abf  fed  fdc 


+  1  fad1  fbc  '’fac’fbd1  ' 


led  to  equations  (4.11),  (4.12)  and  (4.13). 

We  now  prove  a  Schouten-Eisenhart  like  theorem  for  N  =  3 


Theorem  4.3.  is  a  if  and  only  if  for  every  orthogonal 

ennuple  { e  } 


(4.17) 


L  ,  :=  L.  ..  e  1e,  k  =  0 

abc  ljk  a  b  c 


for  distinct  values  of  a,  b,  c  =  1,  2,  3 


Proof.  Let  (e  |  be  an  orthogonal  ennuple  and  define  another 

~cl 

orthogonal  ennuple  { e^ }  by 

a 


e  =  lie  +  ne,  , 

~cl  ~cl  ~u 


(4.18) 


?K  =  -fa  +  M!b  ' 
b 


e  =  e 
~  c  ~c 


for  arbitrary  values  of  the  scalars  ji  and  i>  .  Then  by  (4.17)  we 


(4.19)  L  ...  .  :=  L.  ,  e  1e,: Vk  =  0  , 

abc  ijk  abc 


and  by  using  (4.17)  we  must  have 


Thus,  since  ; i  and  r  ate  arbitrary  it  follows  that 


(4.21) 

L  =  Lkv 

aac  bbc 

(NS) 

for 

distinct 

values  of  a,  b,  c  . 

However,  (4.21)  holds 

for 

a,  b 

t  c  and 

upon  summing  over  a  we 

have 

(4.22) 

>  L  =  2L,  ,  , 

aac  bbc 

a 

(NS 

on  b) 

but 

the  left- 

hand  side  is  identically 

zero 

by  the  doubly  contracted 

differential 

Bianchi  identity.  Thus 

Labc 

=  0  for  distinct  values 

of 

a,  b,  c 

implies  L  =0  for 

aac 

a  t 

c  .  Therefore,  Labc 

=  0 

for 

distinct 

values  of  a,  b,  c  requires 

that  the  space 

be 

conformally  flat.  The  converse  is  obvious  and  our  proof  is  complete. 

The  relative  effectiveness  of  the  Schouten-Eisenhart  and  Ricci- 
Finzi  criteria  can  be  easily  seen.  Suppose  V^tN  >  3)  is  given  and 
it  is  required  to  determine  whethet  is  a  .  By  using  the 

Schouten-Eisenhart  criteria  (Theorem  4.1  or  4.3j  we  would  choose  an 
orthogonal  ennuple  (e^J  and  test  for  =  0  or  =  0  .  If 

R  .  ,  *  0  or  L,  ^0,we  are  done  and  V.,  is  not  a  C„  .  But  if 

abed  abc  N  N 

R  ,  ,  =  0  or  L,  =0  for  je  }  ,  or  even  for  several  different 

abed  abc  ~a 

ennuples,  there  is  no  guarantee  that  V  is  a  since  the  theorem 

requires  that  this  must  hold  for  every  orthogonal  ennuple.  The 
Schouten-Eisenhart  criteria  can  easily  determine  whether 


is  not 


by  a  single  choice  of  ennuple,  but  in  practice  they  are 


inadequate  to  determine  whether  is  a  . 

On  the  other  hand,  the  Ricci-Finzi  criterion  proceeds  by 

attempting  to  solve  the  system  of  partial  differential  equations 

(4.14)  and  (4.15)  for  the  { e  }  .  Upon  converting  these  to  X- 

coefficients  (in  the  notation  used  by  LANDAU/LIFSHITZf 1975] ) 

(4.23)  A  .  .  -  i  ,  , 

abc  abc  acb 

where 


(4.24) 


i  .  «  (A  .  +  A.  -  A  J/2 

abc  abc  bca  cab 


one  can  employ  partial  derivatives,  and  avoid  computing  Christoffel 


symbols  of  VN  . 

Equations  (4.14)  and  (4.15) 

are  equivalent  to 

(4.25) 

A  .  =  0  , 

abc  ' 

(4.26) 

A,,  =  A 

boa  cca 

(NS) 

for  distinct  values  of  a,  b,  c  =  1,  2,  . . . ,  N  .  If  a  solution  of 
these  equations  can  be  found,  then  is  a  .  Moreover,  the 

Ricci-Finzi  criterion  then  produces  the  "nicest"  ennuple  { e^ }  , 

i.e.,  one  which  is  directly  tied  to  the  geometry  of  ,  and  when 

V.  is  a  C,,  a  solution  of  these  equations  is  guaranteed. 

N  N 


h 5.  NORMAL,  GEODESIC,  AND  CANONICAL  CONGRUENCES. 

Definition  5.1.  A  congruence  of  curves  I  in  is  normal  if  the 
unit  tangent  vectors  to  the  curves  of  the  congruence  are  the  unit 
normals  of  a  family  of  surfaces. 


Hence,  if  f  ^  ate  the  components  of  the  unit  tangent  vectors  of 
a  normal  congruence  of  curves,  then  there  exists  a  function 
<t>  :  -*  1R  such  that 

(5.1)  <<>1=‘<)|i=<*>1=+fi- 

In  other  words,  the  unit  tangent  vectors  of  a  normal  congruence 
are  proportional  to  a  gradient. 

Suppose  that  {e  }  is  an  orthogonal  ennuple  in  V,  .  The 
conditions  for  the  congruence  /  defined  by  e,  to  be  normal  are 

<3  -  cl 

well-known  to  be 


(5.2) 


'abc  'acb  -  ^  ' 


where  b  /  c  and  b,  c  t  a  .  See  EISENHART[ 1949 ] ,  p.  115. 


Definition  5.2.  A  geodesic  congruence  of  curves  in  is  one  such 

that  each  of  the  curves  of  the  congruence  is  a  geodesic.  That  is,  if 
f  are  the  components  of  the  unit  tangent  vectors  of  a  geodesic 
congruence,  then 

(5.3)  1.  ^0  =  0  . 

If  (e^)  is  an  orthogonal  ennuple  in  ,  then  the  conditions 


that  e  defines  a  geodesic  congruence  I  ate 


where  b  =  1,  2,  3  .  See  EISENHARTf 1949 ] ,  p.  100.  For  e  to 

-.3 

define  a  congruence  of  curves  that  is  simultaneously  normal  and 
geodesic  it  is  necessary  and  sufficient  that 

(5.5)  e  .  .  =  e  .  .  . 

ai,]  a] , i 

See  EISENHARTf 1949 1 ,  p.  117. 

Suppose  that  £  defines  a  congruence  of  curves  denoted  by 

/  ^  *  .  We  will  construct  two  new  congruences  and  with  unit 

tangent  vectors  A  and  p  ,  respectively,  using  the  method  of 

RICCI [ 1918]  . 

Let  f i  be  the  components  of  f  and  define 


( 5.6 ) 


i-l 


:=  ±  «f; 


1-i 


Consider  the  system  of  equations 


(5.7) 


I'f 

'  (X 


- 1 


i] 


0  , 


'•■gjjK*  +  nfj  =  o  , 


where  and  u  are  scalars  and  C  ate  the  components  of  a 

vector.  The  equations  (5.7)  have  the  determinant 


In  the  following  discussion  it  is  convenient  to  number  the 
congruences  in  a  mannet  which  corresponds  to  oui  labelling  of  the 
ennuple  vectors  (e  }  . 


a 


which  yields  a  second-degree  polynomial  in  <•>  .  The  roots  of  this 
polynomial  are  real  and  when  they  are  inserted  into  (5.7)  determine 
two  real  congruences  of  curves  I  1  and  /  ^  having  unit  tangent 
vectors  \  and  p  ,  respectively.  If  the  roots  are  unique,  then  A 

and  p  are  uniquely  determined.  Regardless  of  whether  A  and  p 

are  unique,  they  are  said  to  be  canonical  with  respect  to  f  ,  and 

I  ^  and  I  1  are  said  to  be  canonical  with  respect  to  I  ^  .  An 
example  of  the  calculation  of  canonical  congruences  is  given  in  the 
appendix . 

We  now  determine  conditions  on  the  rotation  coefficients  of  an 
orthogonal  ennuple  (e^)  when  the  vectors  of  th°  subset  (e^) 

'b  /  a)  are  canonical  with  respect  to  e  .  From  '3.3)  it  is  easy 

to  see  that 


d  i  f  i  g  j 


Applying  the  definition  of  >: 


m 


where  i,^  is  the  root  of  (5.8)  corresponding  to  and  is  the 

scalar  from  (5.7!  corresponding  to  .  Multiplying  (5.10)  by  e  ^ 

with  c  t  a  ,  b  and  contracting  on  3  we  have 


(5.11) 


abc  acb 


=  0  , 


for  b  t  c  and  a  /  b  ,  c  .  Equation  (5.10)  also  implies 


(5.12) 


b  abb  ' 


(NS) 


and 


(5.13) 


-  I  . 

Hb  2  ’baa  ' 


(NS) 


where  a  /  b  . 

In  the  following  we  will  need  to  know  conditions  on  the  rotation 
coefficients  of  the  orthogonal  ennuple  (e  j  when  e  defines  a 

normal  congruence  and  the  vectors  of  the  subset  { e  }  (b  /  a)  are 

canonical  with  respect  to  e  .  Combining  (5.2)  and  (5.11)  we  have, 

in  this  case, 


(5.14) 


abc 


=  0  , 


b  /  c  and  a  *  b  ,  c  .  Rod:  iguer.'  s  formula  McCONNELL  |  1931  ] ,  p.  216 
and  (5.14)  imply  that  the  v-  t  : :  t  a>  are  principal 


directions  on  the  surfaces  n  : :m )  t 

‘  ha 


M  :  <•  ••••«■>  1  .  ih  /  a  )  ate 


the  principal  curvatures 


1  h  /  a 


and 


* — t  |r-4 


inr  v.-r'WV'V.rtf’  'V-  ’V  W  "J"  "J" 

. . . .%vJ 

:$ 


(/^(b  ?  a)  are  the  components  of  the  curvature  vector  of  e^  in 
the  subset  (e^l  (b  /  a)  . 

We  now  have  the  machinery  to  give  our  promised  proof  of  Theorem 
1.2  of  Chapter  III,  the  generalized  Dupin  theorem^ ^  . 

We  prove  only  sufficiency  since  the  necessity  was  proved  in 

Chapter  III.  Choose  two  of  the  families  of  surfaces  and  call  them 
N,  and  .  Let  f.  be  the  components  of  the  unit  normal  f  to 

1^  ,  and  denote  by  I  ^  the  congruence  of  curves  defined  by  f  .  We 

employ  Ricci's  method  to  construct  two  congruences  I  ^  and  r~  with 
unit  tangent  vectors  a  and  p  ,  respectively,  that  are  canonical 

with  respect  to  I  ^  .  Since  '■  and  p  are  principal  directions  on 

'  ,  ,  and  s,  and  v intersect  in  a  line  of  curvature,  then  either 
v  or  p  is  a  principal  direction  on  •  Without  loss  of 

generality  we  may  choose  it  to  be  a  .  Since  and  are 

orthogonal,  c  must  lie  in  the  tangent  plane  to  ,  and  since  f 

is  peipendiculat  to  '  ,  it  is  a  principal  direction  on  > ^  . 

Moieover,  p  is  perpendicular  to  both  a  and  f  ,  and  hence  must  be 

the  unit  normal  to  .  But  since  f  and  a  are  principal 

directions  on  ' ,  they  are  canonical  with  respect  to  p  and  by 

construction  1  and  are  '-anomcal  with  tespect  to  f  .  If  we 


-  .v/ 


A  version  of  th*  fcl ’ owing  matet 
MOORE  ZUND( 1986 ] . 
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will  appear  m  Tensor,  N.S.  See 


s."  v 


>  V* 
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label  our  ennuple  { }  in  the  order  {<\,p,f}  ,  we  may  express  this 
in  terms  of  rotation  coefficients  as 


(5.15) 


•>  -,  ,  =  0  , 
3ab 


where  a  ^  b  and  a,  b  ^  3  ,  and 


(5.16: 


i.  ,  =  0  , 

2ab 


where  a  ^  b  and  a,  b  /  2  .  But  by  the  skew-symmetry  of  the 
rotation  coefficients  in  the  first  two  indices,  (5.15)  and  (5.16) 
imply  that 


7.-0 

abc 


for  distinct  values  of  a,  b,  c  .  However,  (5.1"')  are  necessary  and 
sufficient  for  conditions  that  all  the  congruences  of  curves 
associated  with  the  ennuple  are  normal.  Thus  is  a  normal 

congruence.  Hence,  there  exist  a  third  family  of  surfaces  2^ 
orthogonal  to  both  2^  and  2^  .  This  completes  the  proof. 


Finally,  we  examine  now  canonical  and  normal  congruences  behave 
under  conformal  maps. 


Theorem  5.1.  If  and  are  conformally  related  and  {<\,p,?} 
is  an  orthogonal  ennuple  in  with  A  and  p  canonical  with 
respect  to  f  ,  then  the  conformal  image  of  this  ennuple  in  is 
defined  by 


p 

H 

m 


“  a 

y*y&i 

>>>] 


«  V  /  J< 


v  v 
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>V«VS 
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and  X  ,  p  are  canonical  with  respect  to  f 


Proof.  Let  i  ,  and  .  be  the  rotation  coefficients  of 
abc  abc 

the  respective  ennuples  in  and  .  From  (3.12)  we  have  that 


(5.19) 


2c  r 

^abc  -  e  ^abc 


for  distinct  values  of  a,  b,  c  .  Since  X  and  p  are  canonical 
with  respect  to  §  we  have 


''BIZ  +  '’321  0  • 


(5.20)  7 312  +  7 321  =  0  ' 

and  (5.19)  implies  that 

(5. 2D  7312  +  ,321  =  0  . 

Thus  X  and  p  are  canonical  with  respect  to  f  . 

Theorem  5.2.  Let  and  7^  be  conformally  related.  If  is  a 

normal  congruence  of  curves  in  with  unit  tangent  vectors  f  , 

then  the  conformal  image  is  a  normal  congruence  of  curves  in 

with  unit  normal  f  where  f  is  defined  by 


$ 


i 

A 

/  V  V 


1 


*  1 


Proof.  Use  Ricci's  method  to  construct  congruences  and  I ^ 

with  unit  tangent  vectors  x  and  p  ,  respectively,  which  are 

canonical  with  respect  to  .  Then  we  know  that 


(5.22) 


1 312  ~  "*  321  0  ' 


and  from  (3.12)  we  have  that  t  ,  =  e  i  .  for  a,  b,  c  distinct. 

abc  abc 


(5.23) 


1 312  1 321  0 


f  is  the  unit  tangent  vector  of  a  normal  congruence  I  ^  in 


S6.  A  ROTATION  COEFFICIENT  FORMULATION  OF  THE  CAYLEY-DARBOUX 
EQUATION 


Let  f  define  a  normal  congruence  of  curves  in  V3  .  If  X 
and  p  are  canonical  with  respect  to  f  then 


1  312  *  1 321  "  0  ' 


where  1  ,  are  the  rotation  coefficients  of  (X,p,f)  •  The 
abc  ~  ~  ~ 

condition  that  all  the  congruences  of  an  ennuple  be  normal  is 


(6.2) 


for  distinct  values  of  a,  b,  c  .  The  skew-symmetry  of  the  rotation 
coefficients  in  the  first  two  indices  means  (6.2)  will  hold  when 
(6.1)  holds  and  that 


(6.3) 


1 123  :=  Ai,jp1^  =  ° 


Since  f  defines  a  normal  congruence,  there  exist  <t>  :  ->  (R 


such  that 


(6.4) 


^  ^  • 


If  (6.1)  is  given,  then  (6.3)  implies  that  there  exist  P  , 
9  :  V3  -»  IR  such  that 

I  *i  =  p'\  ’ 

(6.5) 

1  0.  =  Up.  . 
l  l 

Thus  { <#> , vp , 0 }  is  a  triply  orthogonal  system  of  coordinates  in  . 
Hence,  when  (6.1)  is  given,  (6.3)  is  equivalent  to  the  Cayley- 
Darboux  equation  of  Chapter  III,  5  3. 


57.  HOTINE'S  CONJECTURE 

As  stated  in  the  introduction  to  this  Chapter,  Martin  Hotine 
conjectured  and  claimed  to  prove  that  any  function  <*>  :  E^  -*  (R  with 
non-vanishing  gradient  could  be  a  member  of  a  triply  orthogonal 
coordinate  system.  We  quote  from  HOTINEf 1966b] ,  pp.  196-198,  with 
minor  changes  to  agree  with  the  notation  and  numbering  of  this 
dissertation. 


.  W.V.V  v. V.  / 
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"In  HOTINE[ 1966a ] ,  it  was  shown  that  the  gradient  equation  of  a 
scalar  <t>  , 

(7.1)  ♦.  =  •P?i 

(f ^  a  unit  vector),  can  be  transformed  conformally  with  scale  factor 
<#>[=  e°  ]  to  a  Riemannian  space  in  which  the  f  ^  become  tangents  to  a 
family  of  geodesics  and  the  ^-surfaces,  that  is  the  surfaces  over 
which  <t>  is  constant,  become  geodesic  parallels.  It  followed  that 
the  metric  of  the  curved  Riemannian  space  (denoted  by  hats)  can  be 
written  in  the  geodesic  form, 


(7.2)  ds2  =  aapdxad xP  +  d <t>2  (a,p  =  1,2) 

( EISENHART[ 1949 ] ,  p.  57;  WEATHERBURN( 1938 ) ,  p.81).  Using  the  same 
coordinates  (xa,<p)  the  metric  of  the  untransformed  space  can 
accordingly  be  written  as 

(7.3)  ds2  -  a^dx1  dxP  +  e  ^d^2  (a, ft  «  1,2)  . 

"The  components  of  a^  can  of  course  contain  <P  ,  but  for 
different  constant  values  of  <p  will  also  be  the  surface  metrics  of 
the  ^-surfaces.  If  |a|  is  the  determinant  of  the  metric  of  the 
^-surface  passing  through  a  point,  then  it  is  clear  from  (7.3)  that 
the  determinant  of  the  three-dimensional  metric  at  that  point  is 


(7.4) 


Consequently  the  associated  tensor  in  three-dimensions  is 
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■  AAAA  AA  A 


A  AW. 


(7.5) 


(r,s  =  1,2,3) 


r  -  «  -w.  - *  1 
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rs 


(a' 


all 


2n 
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in  which  a”^  is  the  associated  tensor  of  the  surface.  This  can 
easily  be  verified  by  writing  out  the  metric  tensor  in  full. 

"Now  any  curvilinear  coordinate  system  in  three-dimensions 
implies  the  existence  of  three  scalars,  or  coordinates,  whose 
gradient  vectors  are  not  coplanar.  A  coordinate  line  may  be  defined 
as  a  line  along  which  only  one  of  the  scalars  varies,  the  other  two 
being  constant.  Each  coordinate  line  must  accordingly  be 
perpendicular  to  the  gradient  vectors  of  the  other  two  coordinates. 
The  x^-coordinate  line  is  perpendicular  to  the  gradient  of  <p  , 

which  from  (7.1)  is  in  the  direction  fr  normal  to  the  ^-surface, 

1  2 
so  that  the  x  -coordinate  line  (and  similarly,  the  x  -coordinate 

line)  must  lie  in  the  ^-surface,  and  ( x1 , x^ )  can  therefore  be 

considered  as  surface  as  well  as  space  coordinates. 

It  is  apparent  from  the  absence  of  (a, 3)  components  in  the 

metric  (7.3)  that  the  ^-coordinate  line  is  perpendicular  to  both  the 
1  2 

x  -  and  x  -coordinate  lines,  so  that  6  is  the  direction  of  both 

b  r 

the  ^-coordinate  line  and  the  gradient  of  <p  .  Consequently  the 

gradient  vector  of  x^  ,  considered  as  a  scalar  must  lie  in  the 

2 

surface,  because  it  must  be  perpendicular  to  both  the  x  -  and 

2 

^-coordinate  lines.  The  gradient  vector  of  x  must  similarly  lie 

in  the  surface.  We  cannot  yet  say  however  that  the  x^-  and 
2 

x  -coordinate  lines  are  orthogonal  or  that  coordinates  can  be  found 
which  would  make  them  orthogonal  within  the  framework  of  the  space 
metric  (7.3). 


/a 


•  V V  i 
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"We  assume  nevertheless  that  the  coordinates  (x  ,x  )  in  the 
metric  (7.3)  are  definable  as  scalar  functions  of  position  throughout 
some  region  of  flat  space,  in  other  words  that  they  can  be  expressed 
as  some  function  of  Cartesian  coordinates  (x,y,z)  independently  of 
the  definition  of  N  ,  in  accordance  with  the  usual  meaning  of 
coordinates  in  a  Riemannian  flat  space.  In  that  case  we  can  write 


Cp  =  ( x  ) 
*r  r 


where  p^  is  the  unit  vector  in  the  direction  of  the  gradient  of  x1 
and  C  is  the  modulus  of  the  gradient.  Evaluating  C  in  the  metric 
(7.3)  gives 


C2  =  grS!CPt.)(Cps)  =  grS(x1)r(x1)s  =  g11  =  a11 


and  finally 


(7.6) 


c1)  =  vOT  p 

r  a  *r 


But  pr  is  the  unit  normal  to  the  x-surface  passing  through  the 
point  under  consideration,  that  is  the  surface  over  which  x^  , 
considered  simply  as  a  scalar,  is  constant.  Equation  (7.6)  is 
accordingly  in  all  respects  similar  to  (7.1).  By  making  a  conformal 
transformation  to  another  curved  space  with  conformal  factor  a2^  , 
the  p^  will  transform  to  geodesics,  and  exactly  as  in  (2)  we  can 
write  the  metric  of  this  second  curved  space  in  the  geodesic  form 


,2  ,  i  ,6  .,1.2 

ds  =  a  .  dx  dx  +  (dx  ) 

T  r  > 


s*: 


Transforming  back  to  the  original  flat  space,  we  have 

(7.7)  ds2  =  — ^y-  dx'dx'’  +  -jj  ( dx1 ) 2 

a1  a1 

as  the  metric  of  the  flat  space.  We  do  not  yet  know  what  the  other 
"2 

two  coordinates  dx  may  be  in  this  metric.  The  two  metrics  (7.3) 
and  (7.7)  are  however  alternative  ways  of  expressing  the  same  space 
and  one  must  transform  into  the  other. 

"Since  x^  is  a  scalar  function  of  position  in  space,  it  gives 
rise  to  a  family  of  surfaces  whose  metric  can  be  expressed  in  the 
following  forms  putting  dx1  =0  in  (7.3)  and  (7.7) 

,2  ,,2.2  -2o  ,  2 

ds  =  a22(dx  )  +  e  d<f> 


or 


,  2  a->6  ,Ai  ,~f> 

ds  =  — jy-  dx  dx  . 

a1 

These  two  invariant  forms  of  the  line  element  at  a  point  of  the 
surface  hold  not  only  over  one  particular  surface  but  also  over  the 

whole  family,  so  that  the  coordinates  x1  are  either  the  same  as 

2  2 
(x  ,<*>)  or  can  be  transformed  to  (x  ,<t>)  .  We  can  accordingly 

rewrite  (7.7)  as 

i -i  q  \  j  2  2.2  1  ,  ,  1.2  -2n,  2 

(7.8)  ds  =  a22<dx  )  +  -yy  (dx  )  +  e  d<#>  , 

a 
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which  is  triply  orthogonal  in  the  coordinates  (x  ,x  ,<p)  .  Comparing 


this  with  (7.3)  we  have  a^  =  0  ,  which  must  be  so  since  the 
coordinates  are  orthogonal;  and  a^  =  — jy  ,  which  is  so  since 


"We  get  the  same  result  by  considering  the  gradient  of  the  other 
coordinate  . 

"In  aiming  at  this  result,  we  have  merely  assumed  that  the 
1  2 

coordinates  (x  ,x  )  in  the  metric  (3)  are  scalar  functions  of 
position,  without  otherwise  restricting  them  or  the  form  of  <p  . 
Classical  doctrine  on  the  subject  asserts  nevertheless  that  <p  must 
satisfy  a  third-order  partial  differential  equation  known  to 

Eisenhart  and  others  as  the  ( Cayley ) -Da rboux  equation.  If  the  above 
reasoning  is  correct,  then  the  (Cayley) -Da rboux  equation,  which  is 
shown  below  to  be  equivalent  to  one  of  the  six  conditions  of  flat 
space,  must  be  an  identity,  in  which  case  it  expresses  a  relation 
between  and  the  form  of  the  ^-surfaces.  In  the  main  geodetic 

application,  this  would  be  a  hitherto  unsuspected  relation  between 
gravity  and  the  form  of  the  equipotential  surfaces.  It  is 

accordingly  of  considerable  importance  to  resolve  this  question  one 
way  or  the  other. 

"The  remainder  of  this  paper  assumes  the  existence  of  a  triply 
orthogonal  system  derived  from  a  scalar  <p  ,  and  on  that  assumption 
works  out  its  properties,  including  several  which  do  not  seem  to  have 
been  formulated  before,  at  any  rate  in  the  compact  form  now  given. 
If  the  classical  view  is  correct,  these  results  are  valid  if  4> 


satisfies  the  (Cayley) -Da rboux  equation. 


In  the  next  section  we  will  analyze  Hotine's  conjecture 
employing  the  rotation  coefficients  formalism  introduced  in  ^6. 

%8.  CONFORMAL  MAPPING  AND  THE  CAYLEY-DARBOUX  EQUATION 

Antonio  Marussi  was  aware  of  the  weaknesses  of  Hotine's  argument 

4 ) 

and  makes  the  following  comment  in  MARUSSI [ 1985 ] ,  p.  133.  We  quote 
it  using  our  notation. 

”We  now  establish  an  extremely  important  fact.  The  change  in 
curvature  enf..  depends  only  on  the  position  and  on  the  direction  of 
the  normal  to  the  surface,  and  not  on  the  direction  of  the  section  on 
it.  Since  it  is  the  same  for  all  these  directions,  it  follows 
therefore  that  the  directions  of  principal  curvature  are  conserved  in 
the  representation.  Thus,  if  a  family  of  surfaces  is  not  of  Lame's 
type  in  ,  i.e.,  it  does  not  belong  to  a  triply-orthogonal  system, 
then  neither  can  it  transform  in  .  Since  the  family  of 

equipotential  surfaces  of  the  Earth's  gravity  field  is  not  of  Lame 
type  (for  it  to  be  so,  ( Cayley )-Darboux's  third  order  partial 
differential  equation  would  need  to  be  satisfied),  then  neither  can 
be  any  of  its  transforms  in  a  conformal  representation;  there  is  thus 
no  possibility  of  reducing  the  study  of  the  Earth's  potential  field 
into  a  triply  orthogonal  coordinate  system." 

We  will  state  the  above  argument  in  more  precise  mathematical 
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form  and  prove  it. 


Theorem  8.1.  Let  and  be  conformally  related.  Let  I 
denote  a  congruence  of  curves  in  .  If  /  is  not  a  normal 
congruence,  then  its  conformal  image  /'  in  cannot  be  a  normal 
congruence. 

Proof.  The  theorem  follows  by  using  Ricci's  method  to  find  two 
congruences  canonical  with  respect  to  l  and  by  using  equation 
(3.12)  and  the  rotation  coefficient  criteria  in  55  for  a  congruence 
to  be  a  normal  congruence.  See  the  proof  of  Theorem  5.2. 

Theorem  8.1  and  our  derivation  of  the  Cayley-Darboux  equation  in 
Chapter  III,  52  show  chat  Hotine's  conjecture  is  false.  In  511,  we 
will  further  critically  analyze  Hotine's  argument. 

59.  CONDITIONS  FOR  THE  CONFORMAL  IMAGE  OF  A  NORMAL  CONGRUENCE  TO  BE 

A  GEODESIC  NORMAL  CONGRUENCE 

In  this  section  we  examine  conditions  on  the  conformal  function 
a  that  determine  whether  the  conformal  image  of  a  family  of  surfaces 
in  can  be  a  system  of  geodesic  parallels  in  .  (A  system  of 
geodesic  parallels  is  a  family  of  surfaces  such  that  the  unit  normals 
to  the  surfaces  define  a  geodesic  normal  congruence  of  curves.) 

Theorem  9.1.  Let  and  be  conformally  related  and  let  n  be 
the  conformal  function.  If  0  :  -*  IP  is  a  smooth  function  and  we 


define  a  vector  field  f  by 


(9.1)  0 1  =  i f  i 

with  y  >  0  ,  then  f  defines  a  normal  congruence  of  curves,  viz. 
the  curves  normal  to  the  family  of  surfaces  <f>  =  constant.  If  X 
and  p  are  canonical  with  respect  to  f  ,  then  the  necessary  and 
sufficient  conditions  that  the  conformal  image  of  f  ,  f  ,  defines  a 
system  of  geodesic  parallels  are  given  by 


(9.2) 


4<a 


A 


, 

/ 1> 


=  o  , 
=  0  , 


where 


A 


:= 


i 


/a 


:  =  a  .  X 

l 


and 


Proof.  From  (9.1)  we  have 


(9.3) 


<p . 

i .  ]  r  i  i »  d 


but  <t> ^  is  symmetric  in  i  and  j  so 


(9.4) 


*(fi,j  *  ?j,i’  =  Uifj  "  */i 


Multiplication  by  f and  contraction  on  j  y 


)  . 

lelds 


(9.5) 


“i/'lh-Vi,1 


where  i>  = 


i  *  *".  <.  « 


The  conformal  image  of  f  ,  f  ,  is  defined  by 


(9.6)  =  e  fi  , 

and  a  straightforward  computation  shows  that 

l9-7)  r i 3  ‘  - ",  *  fi‘Vf  • 

Combining  (9.7)  and  (9.5)  gives 

(9.8)  >*>f.  =  <#>.  -  /cf  .  -  ^>(n.  -  6  .a  ) 

i ;  j  l  /f  i  i  l  /£ 


If  we 

write  a 

and  *p .  in  terms  of 

l  l 

P  -o  . 
i  /P 

+  ?  •  o  /c 

i  A 

and  =  X.^  +  P^/p 

i 

(9.10) 

^i;/3  =  (*/X  ~ 

A  .  (7  - 

1  /X 


'V  -  <’V'’i 


If  the  equations  (9.2)  hold,  then 


(9.11)  f  •  j  =  0  , 

1 '  J 

and  f  defines  a  geodesic  normal  congruence.  If,  on  the  other  hand, 

f.  #  .(fJ  =  0  ,  then  by  the  linear  independence  of  a  and  p  equations 
i ;  i  ~ 

(9.2)  must  hold.  The  proof  is  complete. 

Corollary  9.2.  Theorem  9.1  remains  true  even  if  and  p  are  not 
canonical  with  respect  to  f  . 


-a.  ’ -  ■  V V_\  V.'.'-.tV.VAts.V.V.VAVulVA  A 


NB:  in  this  corollary  (a, (>,<;}  must  be  orthogonal. 


I.  If  c  =  e‘  ,  then  equations  (9.2)  are  satisfied. 


/  /  V' 

%  v  A 

< .  ■r.-'L 

•vx.** 

/>V. 

£&4 


Theorem  9.4.  The  differential  system  given  by  (9.2)  is  completely 
integrable . 


%•  ’  s. 
^  * 

■M 
'  - 


Proof.  Suppose  we  choose  •*>  =  eT  .  Then  (9.2)  becomes 


(9.12) 


(r  .  -  a.)A  =  0  , 


(t.  -  a.),.  =  0  , 


y- 


and  if  f  =  t  -  a  ,  then  we  have 


(9.13) 


f^1  =  0  , 
f^1  =  0  . 


Let  =  (e  }  a  =  1,  2,  3  .  Equations  <9.13’  may  •  *■■■ 

expressed  in  the  form 


(9.14) 


f.e  :=  f  =  0 
i  a  a 


where  a  =  1,  2  .  These  equa*-  i~r.- 
only  if 


(9.13) 


where  a,  h 
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<9'16»  £/a/b  -  Va  '  £<’cab  '  Wf/c  • 


but  since  f  =  0  when  c  =  1,  2  we  have 


(9.17) 


f/a/b  f/b/a  =  (^3ab  ~  73ba)f/3  ' 


and  the  normality  of  the  congruence  defined  by  f  =  e,  then  requires 


m 


(9.18) 


f/a/b  f /b/a  "  0 


This  completes  the  proof. 


§10.  CONFORMAL  MAPPING  OF  A  V2 

In  this  section  we  consider  conformally  related  two-dimensional 
spaces.  We  do  this  because  much  of  geodesy  is  two-dimensional,  and 
also  to  correct  a  formula  in  MARUSSI ( 1985 ] . 

Let  V2  be  a  two-dimensional  Riemannian  space  with  metric 
tensor  a^  ,  Christoffel  symbols  r^1  ,  curvature  tensor  > 
Ricci  tensor  R^  ,  scalar  curvature  R  ,  and  Gaussian  curvature  K  . 
Greek  indices  will,  of  course,  range  from  1  to  2. 

Let  and  V2  be  conformally  related  with  conformal  function 
o  .  Then  the  metric  tensor  of  V2  is  given  by 


Hi 

w® 


M  *  H  »  »>J 

■v.'  % 


(10.1)  a  :=  e  a  , 

of) 


the  Christoffel  symbols  are 


(10.2)  r  1  :=  f  7  +  f,  V  +  S  V  -  a  ^a^a  , 

up  up  <>  p  p  o  up  o 


(10.5)  R  :=  e  20 (R  +  A2o)  . 

Since  any  is  an  Einstein  space  we  have  that  the  Gaussian 

curvature  of  V2  is 

(10.6)  K  =  e_2°(K  -  &2° )  • 

Suppose  that  A  defines  a  congruence  of  curves  F^  in  V2  . 
If  are  the  components  of  A  ,  then  n  is  given  by 

(10.7)  Pnm,',n\ 

1 

defines  a  congruence  /  2  perpendicular  to  F^  ,  where  tia  is  the 
contravariant  Levi-Civita  dualiEor.  See  McCONNELLt 1931 ) ,  p.  167. 

The  vectors  A  and  p  satisfy  the  Frenet  equations 


McCONNELLt 1931),  p.  185 


We  also  have  a  similar 


where  is  the  geodesic  curvature  of  . 

pair  of  equations: 


(10.9) 


P  * 

p  =  -K  A  , 

'a,p'  go' 

x  P  * 

a,p  g  'o 


where  is  the  geodesic  curvature  of  the  congruence  .  We 

wish  to  exhibit  similar  formulas  for  the  conformal  images  of  and 


r 2  in  V2  ' 


The  conformal  image  of  X  is  defined  by 


(10.10) 


X  =  e  X  , 
a  a 


and  the  conformal  image  of  p  is  given  by 


(10.11) 


a 

=  e  P„ 


By  a  straightforward  computation  we  see  that  the  analogous  formulas 
to  (10.8)  are 


(10.12) 


X  X 

n  ;P 


P 


-o  I  \ 

(>  e  ( x  n  +  o  ) 

"  g  /p 


p„  .nx' 

o  ;P 


-X  e 

(t 


(k  +  a  ) 

g  /p 


where  o 


/p 


Thus  the  geodesic  curvature  of  ('  is  given  by 
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vVvvV 

‘v 


f??7y^ 

(  1  *  »  **  •  -  "1 


"  V 


A>\  > 
u'W’uVk' 


Similarly,  we  have  formulas  analogous  to  (10.9), 


(10.14) 


{; 


P«;PP 


-e"’(Kg*  +  aA)Xn  ' 


.  P  -a  *  * 

A  „p  =  e  (k  +  a  n  )p 

o  ;P  g  A  a 


Hence  we  have  that  the  geodesic  curvature  of  r  ^  is  given  by 


★ 


K 

g 


-a 

e 


"A1 


where  o  .  :=  a  A°  . 

/A  n 

Equation  (10.13)  is  the  correction  to  equation  (4)  of 
MARUSSI[1985]  on  page  150. 


‘ill.  CRITIQUE  OF  HOTINE'S  ARGUMENT 

In  this  section  we  critically  examine  Hotine's  argument  for  his 
conjecture  that  was  quoted  in  ^7. 

In  the  un-numbered  equation  immediately  prior  to  his  (7.6) 
Hotine  has  an  identity  of  the  form 

(11.1)  A1(xi,x-’)  =  g^ 

where  the  (x1}  is  a  coordinate  system  on  (or  E^ )  and  g1-1 

are  the  contravariant  components  of  the  metric  tensor  in  this 
coordinate  system.  However,  (11.1)  is  not  a  general  tensor 
expression,  but  merely  an  algebraic  identity  that  holds  only  at  a 
point.  To  see  this  observe  that  the  left-hand  side  of  (11.1)  are 
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scalars  and  the  right-hand  side  are  the  components  of  a  tensor.  In 
order  to  obtain  the  line  element  (7.7),  Hotine  has  defined  a 
conformal  map  of  -*  by 

(11.2)  9i.  -  l1(x1.x1)9ij  -  9n5i:j  • 

However,  the  second  equality  in  (11.2)  holds  only  at  a  point  and  to 
have  a  valid  definition  of  a  conformal  map  the  conformal  function 
must  be  arbitrary.  Hence  (11.2)  represents  a  non-trivial 
specialization  of  the  conformal  function  a  in 


(11.3) 

Furthermore,  we  have 


2o 


e  g 


ij  * 


(11.4)  glj(k  =  0  , 

where  denotes  covariant  derivative  with  respect  to  the 
Christoffel  symbols  of  g^  ,  and  thus  (11.4)  shows  that  the  g„ 
defined  by  (11.2)  have  the  same  Christoffel  symbols  as  g^  . 

The  remainder  of  Hotine 's  argument  is  entirely  based  on  such  a 
specialization  of  the  conformal  function.  In  the  conformal  geometry 
discussed  in  Chapter  II,  n  is  always  an  arbitrary  function.  Any 
specialization  of  it  generally  makes  it  constant  or  forces  or 
VN  to  be  flat.  Indeed  the  only  known  N  :  3  specialization  of  n 
is  the  concircular  mapping  defined  by 


w§m 

1  V  ^  i 

V  > 
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but  this  is  achieved  only  by  the  introduction  of  a  new  arbitrary 
function  <*>  . 


Hotine  may  have  been  misled  by  two  examples  in  McCONNELL[ 1931 ) , 
which  is  one  of  his  major  references.  The  first  of  these  examples  is 
#4,  p.  197,  which  states: 


If  we  choose  space  coordinates  to  be  orthogonal  and  such 

that  x3  -  0  is  the  given  surface  and  choose  the  u-eurves  on 
the  surface  to  be  the  intersections  of  the  surface  with  the 

x1  -  and  x  -  surfaces,  then  we  have  the  following  relations 
for  all  points  on  the  surface 


(i) 


b 


r 

2 


(ii) 


of  i 


gaft'  9a3 


0;  a 


nf  I  nfi 


g"3  -  o,  g33  . 


'33 


( iii ) 


). 


The  only  thing  explicit  in  McConnell's  statement  is  that  the  xr 
are  orthogonal  curvilinear  (not  Cartesian)  coordinates  for  E^  .  A 
similar  ambiguity  occurs  in  EISENHARTf 1947 ) ,  p.  159. 

Presumably  (i)  is  a  formal  specialization  of 


(11.6) 


,  r 

1 1 

s 


which  is  a  familiar  result,  e.g.,  EISENHARTf 1949 ] ,  p.  2,  that  occurs 
in  the  classical  treatment  of  tensor  transformation  laws.  However, 


vVv 

m 


9? 


(11.6)  is  valid  only  at  a  point  —  a  fact  not  stressed  by 
McCONNELL[ 1931 ]  but  which  is  always  implicit  in  the  tensor 
transformation  laws. 

If  we  consider  (i)  as  a  system  of  differential  equations  we  have 


(11.7) 


1  1  1 

X  =  u  +  c 

2  2  2 

X  =  u  +  c 


3  3 

X  =  c 


12  3 

where  c  ,  c  ,  and  c  are  constants.  Thus  if  the  differential 
equations  hold  at  all  points  of  the  surface,  then  the  surface  would 
be  flat.  This  is  easy  to  see  if  one  examines  the  Lame  equations  in 
Chapter  III,  S3  and  the  usual  formula  for  Gaussian  curvature  when  the 
surface  metric  is  orthogonal,  viz. 


1 

2>/a 


((  a22|l)jl  +  ( 


‘nU’  |i] 


It  can  be  arranged  that  (i)  and  hence  (ii)  and  (iii)  hold  at  a  point 
of  the  surface.  Hence  this  example  is  incorrectly  stated  —  the 
quoted  results  are  valid  only  at  a  point,  not  on  a  surface. 

McConnell's  second  example  is  41,  page  188,  which  is  similar  to 
equation  (11.1)  and  states  that 


(11.8)  A  j  ( u"  ,i/! )  =  a^1 

where  a"^  are  the  contravar lant  components  of  the  surface  metric 
and  u"  are  coordinates  on  the  surface.  The  second  example 


essentially  involves  the  first  since  by  definition 


,  o  d.  '»<*>  (Iu°  <lun  «P 

11.9)  A1(u  ,u  )  »  a  — - - ^  =  n  b^£>b  =  a  . 

1  c)u'1  du  ’ 


Presumably  McConnell's  second  example  was  suggested  by  the 
discussion  of  BIANCHI[1910] ,  pp.  67-69,  which  is  correct.  Since  this 
discussion  is  very  instructive  we  will  briefly  outline  it.  Bianchi 
seeks  to  construct  a  non-singular  change  of  variable 


......  O  — O  (T  ,  1  2  . 

(11.10)  u  — >  u  =  v  (u  ,u  ) 

which  reduces  the  first  fundamental  form  ds^  =  a^du  Mu'3  on  a 

surface  ,  to  a^dv  W*  having 

(H.H)  an  -  i22  ,  S12  -  0  . 

Then  by  (11.10),  the  correct  statement  of  (11.8)  is 


(11.12)  A1(v°,v/?)  =  a*P  . 

This  follows  by  the  definition  of  A^  (see  Chapter  II,  §3) 


(11.13) 


pa  civ 

i  - 

R. 

It 

o 

pn  c7u 

a  - 

dP 

duP 

du 

dup 

du 

which  is  the  classical  tensor  transformation  law.  However,  (11.13) 

holds  not  only  at  a  point  but  in  a  coordinate  chart  of  the  point. 

1  2 

Writing  f  :=  v  and  <i  :=  v  ,  (11.12)  becomes 


(11.14) 


.  <-  -11 
Ajf  -  a  , 


t  ,  -12 

A^((f  ,'l)  =  3  , 


A  jU 


-22 

a 


But  since 
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(11.22)  Ay/  =  0  . 

The  systems  (11.19)  and  (11.21)  are  called  Beltrami  systems  and  they 

are  the  Cauchy-Riemann  equations  on  the  curved  surface  V-  . 

Likewise  (11.20)  and  (11.22)  are  Laplace's  equations  on  V~  • 

Thus  the  required  reduction  of  a^^dxfdxP  to 

2  2 

,\(f,g)  {df  +  dr/  )  is  equivalent  to  determining  solutions  of  a  pair 

of  Beltrami  systems.  This  problem  is  now  other  than  the  construction 

of  an  isothermal  coordinate  system  on  ,  and  in  effect  explains 

our  comment  on  page  26  that  any  is  conformal  to  any  other  . 

There  are  an  infinite  number  of  systems  of  isothermal  coordinates, 

each  system  corresponding  to  an  analytic  function  of  the  complex 

variable  u  +  iu  More  precisely  stated,  if  the  coordinate  lines 
1  2 

u  =  constant,  u  =  constant  then  all  other  isothermal  systems  are 
given  by  the  equations 

1  2 

Re{f(u  +  iu  )]  =  constant, 

(11.23) 

Im(f(u^  +  iu^"))  =  constant, 

1  2 

where  f  is  an  analytic  function  of  u  +  iu  Thus  corresponding 
to  each  isothermal  system  there  is  a  conformal  mapping 


All  of  Bianchi's  work  is  two-dimensional;  however,  Hotine's 
argument  requires  a  more  complicated  construction  in  three 
dimensions,  and  his  proof  involves  a  serious  omission.  In  effect  he 
specializes  the  conformal  function 


(11.24) 


to  obtain  the  line  element  (7.3).  This  is  employed  to  map  the  system 
of  surfaces  <p  «  constant  into  a  system  of  geodesic  parallel 
surfaces.  Since  the  function  <t>  is  not  truly  arbitrary,  it  remains 
to  be  shown  that  such  a  specialization  is  valid.  In  order  to 
complete  Hotine’s  argument,  it  would  be  necessary  to  prove  that  his 
specialization  was  admissible.  In  view  of  the  previously  indicated 
error  there  seems  little  reason  to  attempt  such  a  proof.  Moreover, 
it  seems  highly  unlikely  that  such  a  proof  could  be  done  since  the 
requirement  that  a  map  a  system  of  geodesically  parallel  surfaces 
into  a  system  of  geodescially  parallel  surfaces  forces  a  to  reduce 
to  a  constant,  i.e.,  a  homothety.  Moreover,  Bianchi's  reduction  of 
a^du'W*  to  isothermal  form  was  possible  only  by  the  existence  of 
harmonic  conjugate  functions,  and  is  a  direct  consequence  of  the 
plentiful  supply  of  analytic  functions  of  a  complex  variable.  In 


Hotine's  case  he  requires  a  more  complicated  specialization  of 
g.  .dx1dx-'  to  the  form  (7.3).  However,  in  three  dimensions  the 
situation  is  significantly  different.  The  three-dimensional  analogue 
of  the  Beltrami  systems  (see  HEDRICK-INGOLDf 1925 ] )  does  not  reduce  to 
Cauchy-Riemann-like  systems  and  the  supply  of  analytic  functions  — 
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hence  conformal  maps  —  is  quite  limited.  This,  of  course,  was 
predicted  by  the  Liouville  Theorem  in  Chapter  III,  $3.  For  this 
reason  it  seems  unlikely  that  Hotine's  argument  could  be  completed. 

Finally  we  can  also  consider  Hotine's  conjecture  strictly  from 
the  viewpoint  of  partial  differential  equations,  let  f  :  -»  IR  be 
a  given  function  with  non-vanishing  gradient.  From  partial 
differential  equations  we  know  that  the  system 

g1^ ^  :=  a  1  (f ,n)  =  0 

(11.25)  | 

9ljfiCj  :=  A1(f ,C )  =  0 

has  solutions  n»C  :  -*  IR  .  If  in  addition  we  require  that 

(11.26)  gljC j  -  A1('hC)  =  0  , 

that  is,  define  a  triply  orthogonal  system,  then  the  system 

of  partial  differential  equations  given  by  (11.25)  and  (11.26)  is 
over-  determined  and  need  not  have  a  solution.  Another  equation  on 
f  is  needed  to  ensure  a  solution  and  this  new  equation  is  precisely 
the  Cayley-Darboux  equation  of  Chapter  III.  In  fact,  the 

Cayley-Darboux  equation  is  derived  from  a  system  of  partial 
differential  equations  analogous  to  (11.25)  and  (11.26). 

It  is  now  easy  to  see  the  connection  between  the  generalized 
Dupin  theorem  and  the  Cayley-Darboux  equation.  If  we  have  a  system 
of  surfaces  defined  by  f  =  constant,  then,  it  is  alv/ays  possible  to 
find  another  system  orthogonal  to  it.  However,  the  generalized  Dupin 
theorem  states  that  the  existence  of  a  third  system  orthogonal  to 


both  of  the  others  requires  the  original  two  systems  to  intersect  in 
lines  of  curvature.  Therefore,  we  may  restate  the  Cayley-Darboux 
theorem  in  the  following  manner. 

If  f  :  -»  K  defines  a  system  of  surfaces,  then  in  order  for 

there  to  exist  another  system  of  surfaces  orthogonal  to  that  defined 
by  f  and  such  that  the  two  systems  intersect  in  lines  of  curvature, 
f  must  satisfy  the  Cayley-Darboux  equation. 

In  conclusion  we  have  indicated  three  major  reasons  why  Hotine's 
argument  is  seriously  flawed: 

(i)  it  involves  a  choice  of  coordinates,  and  simplification  of 
original  metric  which  is  valid  only  at  a  point; 

(ii)  it  employs  a  specialization  of  the  conformal  mapping 
function,  but  does  not  establish  the  admissibility  of  this 
specialization; 

(iii)  it  does  not  verify  that  the  geoidal  surface  given  by 
<t>  =  constant  satisfies  the  Cayley-Darboux  equation. 

Any  one  of  these  reasons  would  be  non-trivial  to  rectify,  and  taken 
together  we  feel  they  show  that  Hotine's  argument  is  fatally  flawed. 
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Appendix  -  EXAMPLE  OF  THE  CALCULATION  OF  CANONICAL  CONGRUENCES 

Let  =  IR^  =  E^  and  set 

,  ,.123,  ,1,2  ,2,2 

(A.l)  f(x  ,x  ,x  )  =  (x  )  +  (x  ) 

2 

and  for  the  purpose  of  convenience,  rewrite  (A.l)  as  f  »  x^  +  x 
Then 


(A. 2) 


vf  =  2(x1<x2,0) 


hence  the  unit  normal  to  the  surfaces  f  =  constant  is 


(A. 3) 


e  = 


2f 


1/2  : 


vf  , 


and  (5.8)  becomes 


(A4) 


2f 

1 


3/2 


2f 


3/2  xlx2 


2f 


1/2 


2  f 


3/2  X1X2 


2f 3  2 


2f 1  2 


'1 


2f 

x 


1/2 

2 


2f 


12 

0 

0 


Evaluation  of  this  determinant  yield: 


(A. 5) 


r  1 

,''1  4  '  nA  2 


=  0 


at 


which  implies  that 


Using  the  second  equation  of  (5.7)  the  vector  corresponding  t 
i.>.  =*  0  is  given  by 


(A. 7) 


X  .a\\]  =  0 
1] 


If  we  make  the  change  of  variables 

(A. 9)  x.  =  i  cosd  ,  x,  =  i  sinh  ,  x  ■  x 

l  i.  i 

we  obtain 


(A. 10) 


A  *= 


«o,o,n  , 


(A. 11) 


n  =  (-  sinO,cosO,0)  , 


(A. 12) 


f  =  ( cost) ,  sinO ,  0  ) 


which  ate  the  well-known  basis  vectors  for  cylindrical  coordinates. 
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SUPPLEMENT  (by  J.D.  ZUND) 


$1 .  Introduction 

The  work  in  Dr.  Moore's  dissertation  strongly  suggests  that  Hotine's 
conjecture  is  highly  implausible.  It  shows  that  crucial  steps  in  his 
arguments  are  fatally  flawed,  but  it  does  not  prove  that  the  conjecture  is 
false.  Without  using  Hotine's  equations  we  cannot  obtain  his  results,  and  we 
were  unable  to  derive  his  equations  by  other  means.  Our  goal  was  then  to  find 
a  decisive  part  of  his  argument  which  was  wrong.  We  succeeded  in  doing  this 
about  a  week  after  the  dissertation  was  submitted  to  the  examining  committee 
and  the  Graduate  School  of  New  Mexico  State  University.  In  this  supplement  I 
will  describe  this  conclusive  error  in  more  detail  than  was  possible  in  our 
joint  paper  "Hotine's  Conjecture  and  Differential  Geodesy,"  which  will  be 
published  in  Bui letin  Goodrsique . 

As  mentioned  in  our  introduction  to  this  research  report,  the  major 
oversight  on  Hotine's  part  was  his  failure  to  employ  the  formalism  of  Ricci 
rotation  coefficients.  In  Section  ‘<2  we  will  translate  Hotine's  Cayley- 
Darboux  equation  into  the  language  of  rotation  coefficients.  In  Section  N 3  it 
will  be  shown  that  his  Cayley-Darboux  equation  is  not  equivalent  to  the  true 
Cayley-Darboux  equation  and  moreover  neither  of  these  equations  is  an 
identity.  Finally,  in  Section  *»4  we  will  discuss  the  physical  ramifications 
of  this  result  and  why  MARUSSI  [1985]  (page  133)  was  correct  in  his  1967 
statement  that  "there  is  no  possibility  of  reducing  the  study  of  the  earth's 
potential  field  onto  a  tr lply-orthogonal  system." 

'•2.  Hotine's  Cayley-Darboux  Equation 

In  this  section  we  will  reformulate  H  tine’s  Cayley-Darboux  equation  in 
tne  language  of  rotation  coefficients.  Bef.--.te  do  inn  this,  a  few  general 
comments  on  notation  are  required.  Fust,  >  ut  paper  "Hotine's  Conjecture  and 


Differential  Geodesy"  for  Bulletin  Geodesique  was  essentially  written  in 
Hotine's  notation  which  differs  slightly  from  that  employed  in  Dr.  Moore's 
dissertation.  The  main  differences  are  that  Hotine  did  not  employ  a  comma  to 
denote  covariant  derivatives,  and  we  indicate  partial  derivatives  by  a  stroke 
"I".  Moreover,  Hotine's  conformal  function  m  appears  in  the  dissertation 
as  e“  ,  and  the  triad  (A^r)  in  our  paper  is  denoted  in  the  dissertation  by 
Since  this  supplement  is  intended  to  be  read  in  conjunction  with  our 
paper,  it  will  be  written  in  Hotine’s  notation.  Additional  differences  in 
notation  will  be  noted  where  necessary. 

Second,  references  to  equations  appearing  in  the  dissertation  will 
include  a  page  citation.  References  to  our  Bulletin  Geodesique  paper  are  to 
equations  appearing  in  its  appendix  and  will  include  the  prefix  "A"  in  the 
equation  number.  All  numbered  equations  in  this  supplement  will  include  the 
prefix  "S”  in  the  equation  number.  Finally,  references  in  the  supplement  will 
always  refer  to  those  listed  in  the  bibliography  of  the  dissertation. 

In  HOTINE  [1966b]  (page  202),  and  in  HOTINE  [1969)  (page  114),  Hotine 
announced  two  forms  of  the  Cayley-Darboux  equation:  a  sin  face  equal  ion 


(£)„/' m'  =0  (o,  I1  =  1,  2) 


and  a  s/*m  e  aqua f i on 


(i>  V„S  -  0 

n  rs 


( r,  s  =  1,  2,  3)  . 


In  both  these  equations  the  vectors  a",  iJ1  and  A r ,  |iS  respectively  are 
principal  directions.  The  factor  n  appears  in  the  basic  gradient  equation: 


( S. 3)  N  =  nr  , 

r  r 

and  in  this  equation  N  is  to  be  identified  with  the  geopotential  of  the 
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rotating  Earth;  each  N-surface  N(x  ,x  ,  x  )  =  constant  is  an  equipotential 
surface;  and  *>r  is  the  unit  normal  to  the  N-surfaces.  Hence  physically  n 
corresponds  to  the  magnitude  of  the  gravitational  acceleration.  In  the 

dissertation,  (S.3)  is  written  with  N,  n,  i>  replaced  by  <f> ,  -p,  (see 

(7.1)  page  75) . 

Strictly  speaking,  Hotine's  derivation  of  both  (S.l)  and  (S.2)  are  rather 
dubious.  In  HOTINE  [1969)  (pages  113-114)  few  details  are  given  and  we 
suspect  that  the  missing  steps  include  use  of  equations  which  are  valid  only 
at  a  point  (see  the  dissertation  pages  90-91).  However,  the  real  issue  is  not 
how  he  derived  these  equations,  but  the  fact  that  he  regarded  them  —  in 
particular  (S.2)  —  as  the  genuine  Cayley-Darboux  equations. 

In  order  to  translate  (S.2)  into  rotation  coefficient  language,  we  take 

the  unit  vector  r  to  be  the  third  vector  of  our  triad  { X ^ r } ,  with  the  unit 

vectors  and  p  being  the  first  two  vectors  of  the  triad.  These 

vectors,  or  equivalently  the  congruences  J ^  and  l  ^  having  them  as  tangent 

vectors,  are  assumed  to  be  canonical  with  respect  to  ;>r.  Hotine  did  not  use 

this  terminology,  but  this  important  property  is  implicit  in  his  construction. 

This  notion  is  fully  discussed  on  pages  65-68  of  the  dissertation.  The 

resulting  triad  is  denoted  as  (a  r)  and  explicitly  it  is 

a 

,,  r  r  r. 

We  now  proceed  to  translate  (S.3)  into  rotation  coefficients.  First  we 
covanantly  differentiate  it  to  obtain 


( S .  4  ) 


and  similarly 


N  -  n  r  *  nr 
rs  s  r  rs 


N  =  n  i’  +  m 
sr  r  s  st 


But  since  Nf  is  a  gradient,  both  these  expressions  are  equal,  viz 
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Using  (S.4)  and  (S.5)  we  obtain 


N  =  N  ! 
rs  sr 


since  i> 
obtain 

( S.6 ) 


..  r  s  ,  .  r  s 

N  i>  u  =  (n  i)  +  nn  )  i>  /j  =  n  , 
rs  s  r  rs  /2 


Nsr"r,jS  =  (rY’s  +  ni'sr )|,r,,S  =  m323  ' 

r  r 

is  a  unit  vector,  i.e.,  i*  n  =1,  we  have  n  i>  =  0. 

[  lb 


n/2  ni 323 


Thus  we 


where  n^  =  n^i" 


and  "/"  is  defined  as  in  (A. 2),  or  in  the  dissertation  on 


pages  54  and  58 
yields  the  corresponding  equation 

(S.7) 


r  s  r  s 

Repeating  the  same  procedure  for  N  i>  ^  and  Nsru  ^ 


n/l  ”  m  313 


where  n  =  n  Xs. 

A  s 


_  1 


It  is  convenient  now  to  temporarily  write  '!'  =  ->  so  (S.2)  becomes 


( S . 8 )  *  Am  =  0  . 

rs 

To  convert  this  into  rotation  coefficients,  we  compute  the  derivatives 


4-  =  -  ^  n  , 

r  2  r 
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n/l/2  =  (nrsA  "  +  V  s'J 


However,  we  have 


.  r  s  r  r 

A  s''  =  1122"  +  1 1321’ 


and  since,  by  hypothesis  I  ^  is  normal,  ^  321  f'see  Pa9e  66)  and 

canonical  ^  =  0  (see  (5.11)  page  69),  this  reduces  to  merely 


Thus,  we  obtain 


i  .e. , 


(S.10) 


.  r  s  r 

A  s"  =  ^122"  • 


n  —  n  A  1  j  "^*"^1  -,->n  , 

/1/2  rs  122  /2 


n  A  u  _  n  ..  1 .  „n 

rs  1  /1/2  122  /2 


Now  using  (S.7)  and  differentiating,  we  see  that 


n/l/2  "  1 313n/2  +  m  313/2 


and  by  (S.6)  this  yields 


n/l/2  =  n  "*  3231  313  +  1  313/2'  ' 
Thus,  combining  these  with  (S.10)  we  have  finally 


nrsA  11  ~  n(',3231313  +  1  313/2  ^  323'’  122 }  ' 

If  this  is  substituted  into  (S.9)  we  obtain 


n  ('T3231  313  +  ^  313/2  ',3231  122  2~’323'’313)  0 


which  simplifies  to 


(S.ll) 


1 313/2  1 3231 313  ’ 3231 122  ~  ° 


and  this  is  Hotine's  Cayley-Darboux  equation. 


1 


5  3 .  The  Cayley-Darboux  Equation 

We  now  establish  that  Hotine's  Cayley-Darboux  equation  is  not  equivalent 
to  the  true  Cayley-Darboux  equation,  and  that  neither  this  equation,  nor 
( S . 11 ) ,  is  an  identity.  This  conclusively  refutes  the  conjecture. 

The  first  part  is  easy.  In  the  case  of  a  triply-orthogonal  system  of 
surfaces,  the  respective  surfaces  pairwise  interest  in  the  normal  canonical 
congruences  f^,  f’3  as  described  in  the  triad  scheme  of  52.  Thus,  when 
^312  =  ^321  =  ®  (the  conc^ition  f°r  the  normal  congruences  and  f  ^  to  be 
canonical  with  respect  to  r,),  as  shown  in  the  dissertation  on  page  74,  the 
Cayley-Darboux  equation  assumes  the  simple  form 

(S.12)  ^123  =  0  . 

It  is  now  obvious  that  this  does  not  have  the  same  structure  as  Hotine's 
equation  (S.ll). 

Hotine  claimed  that  his  equation  was  an  identity  by  virtue  of  the  Lame 
equations  (page  25)  which  express  the  flatness  of  Euclidean  3-space.  In  terms 
of  rotation  coefficients  these  are  expressed  by  the  six  equations 


( S . 13 ) 


R  ,  ,  =  0  , 

abed 


1212 

=  o, 

R1213  °' 

R1 223  0  ' 

1313 

=  0, 

R1323  =  °' 

R2 323  =  0  ' 
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By  an  inspection  of  (S.14)  it  is  clear  that  only  one  of  the  Lame  equations  has 
a  form  similar  to  (S.ll).  This  is  ^323'  by  skew-symmetries  is  R3i 32 ' 


R3132  7 313/2  7 312/3 

+  ^ ( 7  f 327  f 13  -  7  f  3 37  fl2 


+  Z7  31f ( 7  f 32  7f23* 
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Expanding  this  expression  using  the  normal  and  canonical  condition 
7 312  =  7 321  =  we  obtain 

(S,15)  0  "  7 313/2  ”  7 3137 323  "  7 323^ 122  "  7123<7  311  ~  7 322 1 

where  the  underlined  terms  are  precisely  those  appearing  in  (S.ll).  This 
shows  that  Hotine's  equation  is  merely  u  piece  of  «  lame  equation,  and  (S.ll) 
does  not  imphj  that  y  ±23  =  ®  which  is  the  true  Cayley-Darboux  equation. 
Moreover,  (S.15)  shows  that  ^23  nee^  n°t  zeco-  hence  the  Cayley-Darboux 

equation  is  certainly  not  an  identity! 


*?4.  Physical  Consequences 

In  (i3  we  showed  that  neither  Hotine's  Cayley-Darboux  equation  or  the  true 
Cayley-Darboux  equation  is  an  identity,  i.e.,  a  consequence  of  the  Lame 
equations  for  the  flatness  of  .  Thus,  in  order  for  to  admit  a  triply- 
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orthogonal  system  of  coordinates  the  true  Cayley-Darboux  must  be  imposed  as  an 
additional  condition.  This  is  presumably  the  reason  for  the  cryptic  comment 
of  Marussi  quoted  at  the  end  of  the  introduction  to  this  supplement. 

The  reasons  for  such  a  conclusion  are  now  obvious,  from  the  complicated 
nature  of  the  Cayley-Darboux  equation  (recall  n6  of  Chapter  IV,  pages  73-74). 
One  can  argue  both  on  mathematical  and  physical  grounds,  and  in  fact  such 
reasons  are  not  truly  independent. 

Suppose  one  considers  the  simple  case  of  the  Earth  rotating  with  a 
uniform  angular  velocity  «.»  and  having  N  as  its  geopotential  function. 
Then  one  has  the  Newtonian  equation 

( S.  16 )  AN  =  -2(/ 

(where  A  is  the  3-dimensional  Laplacian),  and  if  the  function  N  is  to 
define  a  N-surface  of  a  triply-orthogonal  system  of  surfaces,  then  it  must 
also  satisfy  the  Cayley-Darboux  equation 

( S . 17 )  KN  =  0  . 


Equation  (S.17)  is  a  convenient  way  of  writing  the  Cayley-Darboux  equation  (it 
was  not  used  in  Chapter  III,  or  ZUND/MOORE  (1986),  but  it  was  introduced  in 
our  Bulletin  Geodesique  paper).  In  effect,  K  might  be  called  the  Cayley- 
Darboux  operator  which  as  we  know  is  a  complicated  non-linear  third  order  non¬ 
linear  partial  differential  operator.  A  tensor  expression  for  (S.17)  was 
given  in  ZUND/MOORE  [1986],  i.e., 


Mathematically  for  a  triply-orthu  i  :  •.  -.yet era  one  is  faced  with  solving 
the  combined  system  of  equations  (S.16)  and  < S . 1 7 ) .  This  system  is  rather 
unspeakable,  and  there  is  no  reason  to  suspect  that  these  two  equations  are 
even  consistent.  Even  worse  are  the  technical  questions  of  the  existence  and 
uniqueness  of  such  a  combined  linear  and  non-linear  set  of  equations. 
Virtually  nothing  is  known  in  this  respect!  From  the  underlying  physical 
situation,  one  would  like  to  regard  (S.16)  as  the  basic  equation  for  N  with 
(S.17)  being  some  kind  of  geometric  constraint  on  the  function  N.  However, 
mathematically  these  roles  might  be  interchanged!  Even  if  all  these  questions 
are  successfully  resolved  and  (S.17)  is  regarded  as  a  geometric  constraint  on 
the  problem,  there  is  no  guarantee  that  it  will  not  exclude  all  the 
interesting  solutions  of  (S.16).  Thus,  mathematically  the  combined  system  is 
a  horrendous  problem. 

Physically  the  situation  is  much  simpler.  It  is  clear  that  (S.16)  is  the 
significant  equation  and  that  it  is  sufficient  to  determine  the  Newtonian 
geopotential  function  N.  The  theory  is  complete  and  requires  only  him- 
c(/i Kit  inn,  (S.16),  and  there  is  really  no  need  for  (S.17)  in  the  physics  of  the 
Earth's  gravity  field.  Barring  an  unlikely  physical  interpretation  for  (S.17) 
—  which  would  make  potential  theory  into  a  non-linear  theory  (a  non-Newtonian 
theory)  —  this  equation  is  nonsense.  Its  only  purpose  is  to  produce  a 
triply-orthogonal  system  of  surfaces  and  a  'nice'  coordinate  system.  However, 
its  complicated  nature  --  both  mathemat ically  and  physically  —  since  it  could 
exclude  physically  meaningful  solutions  of  (S.16)  —  suggests  that  such 
coordinate  systems  need  not  exist.  We  have  no  doubt  that  if  Hotine  had  known 
that  the  correct  Cayley-Dar  boux  was  :v't  in  identity,  he  would  have  not 
proposed  his  con]ertute.  Stated  siif.--.net  ly,  i.S.17)  is  not  a  physical 
equation,  and  it  does  not  fit  into  the  ti  itui<-  of  Newtonian  potential  theory 
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